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ABSTRACT 

We present the first single-burst stellar population models in the infrared wavelength 
range between 2.5 and 5 fim which are exclusively based on empirical stellar spectra. 
Our models take as input 180 spectra from the stellar IRTF library. Our final single¬ 
burst stellar population models are calculated based on two different sets of isochrones 
and various types of initial mass functions of different slopes, ages larger than 1 Gyr 
and metallicities between [Fe/H] = —0.70 and [Fe/H] = 0.26. They are made available 
online to the scientific community on the MILES webpage. We analyse the behaviour 
of the Spitzer [3.6] —[4.5] colour calculated from our single stellar population models 
and find only slight dependencies on both metallicity and age. When comparing to 
the colours of observed early-type galaxies, we find a good agreement for older, more 
massive galaxies that resemble a single-burst population. Younger, less massive and 
more metal-poor galaxies show redder colours with respect to our models. This mis¬ 
match can be explained by a more extended star formation history of these galaxies 
which includes a metal-poor or/and young population. Moreover, the colours derived 
from our models agree very well with most other models available in this wavelength 
range. We confirm that the mass-to-light ratio determined in the Spitzer [3.6] /rm band 
changes much less as a function of both age and metallicity than in the optical bands. 

Key words: infrared: stars - stars: fundamental parameters - infrared: galaxies - 
galaxies: stellar content - galaxies: structures 


1 INTRODUCTION 

Most galaxies cannot be resolved photometrically into 
single stars due to their large distances. In order to study 
their stellar content, observations of their integrated light 
can be analysed by means of so-called single-burst stellar 
population (SSP) synthesis models which re produce a pop- 
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Conrov & van Dokkum 2012f). 


These models are constructed by populating stellar evolu¬ 
tionary isochrones with stellar spectra following a particular 


inate 


3 from large stellar libraries (e.g.. 
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1: IPrueniel & SoubiranI 200ll: Le Borene et al. 
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200911 . These are compilations of either 


observed or theoretically constructed stellar spectra, ideally 
covering a large range of the stellar atmospheric parame¬ 
ters effective temperature Tefr, surface gravity log(y) and 
metallicity [Fe/H]. Whereas a number of empirical stellar 


libraries exist in the optical wavelength range (e.g.. 

Pickled 

19981: Prugniel & SoubiranI 200ll: Le Borgne et al 

1200 . 3 I: 

Valdes et al.l l2004l: ISanchez-Blazauez et al.l l200fil. Pickles. 


ELODIE, STELIB, Indo-US, MILES), only very few 
empirical libraries are available in the ne ar-infrared (NIR) 
and infrared (IR) wavelength range (e.g., iLancon fc WoodI 
I 2 OOOII . Throughout this paper, we will refer to the wave¬ 
length range until the end of the K-band as NIR, and 
as IR to the range between 2.5 and 5 /rm, in accordance 
with the generally used nomenclature. Mo st SSP models 
employ theor e tical s tellar libraries (e.g., jKuruc i I 1992 I: 
iLeieune et al.l Il997l . 1 19981 : IWestera et al.i I 2 OO 2 I . BaSeL) 
in this wavelength range. These are, however, generally 
hampered by low resolution and uncertain prescriptions 
of theoretical model atmospheres. On the other hand, the 
caveats of empirical stellar libraries are their natural limits 
in parameter coverage. Since most of the stars of the stellar 
libraries are observed in the Solar neighbourhood, their 
metallicity coverage is rather small around the Solar value. 
However, the quality of the resulting SSP models depends 
very much on the parameter coverage of the input stellar 
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libraries. The NASA Infrared Telesc ope Facility ( IRTF) 
spectral library dCushing et al.l l2005l : iRavner et al.i l2009l ~) 
is the only empirical library available at the moment that 
extends towards the IR with a reasonably good stellar 
parameter coverage (see Section 

There are various ways in which observed galaxies can be 
studied with the help of SSP models. One possibility is to 
compare their observed colours with models. Another option 
is to focus on absorption line strength indices visible in both 
the recorded and the modelled spectra fe.g.. IWorthev et al.l 
19921: IVazdekis et al.lll99l l2003l: iThomas et al.ll200l l2005l: 


Cervantes fc Vazdekisi ' 20091 : iThomas et al.l I 2 OI 0 II or even 
carrying out a full spec t rum fitting approa ch (see e.g. 
ICid Fernandes et al.l I 2 OO 5 I : iKoleva et al.l I 2 OO 9 I I . In the lat¬ 
ter case, the whole spectral energy distribution of an ob¬ 
served unresolved stellar population is reproduced by a best¬ 
fitting combination of model spectra. This way, valuable con¬ 
straints on the underlying stellar populations and their star 
formation histories (SFH) can be drawn. 

The general lack of observations of both stars and galaxies 
in the NIR and IR wavelength range gets also reflected by 
the small number of models available so far in this regime. 
Models like the Galaxy Ev olutionary (GALEV) synthesis 
models (iKotulla et al.ll2009li. the Flexible Stella r Population 
S ynthesis (FSPS) models ijConrov et al.l 120091'). the m odels 
of iBruzual fc Chariot! ll20oji . the models of MarastonI (l2005l l 
and others which are based on the theoretical B aSeL-library 
llLeieune et al.lfT997l . Il998l : IWestera et al.l[2002h are the only 
on es providing model spectra up to the IR. The latter ones 
by iBruzual fc Chariot! (l2003l l a re also available based on 
the empirical STELIB library llLe Borgne et all l2003l l for 
the optical wavelength range. A number of other models 
exist in the NIR wavel ength range up to t h e end of the 
K band. For example, iMouhcine fc LanconI ll2002li calcu¬ 
lated m odels in this regime base d on the empirical stellar li¬ 
brary of iLancon fc WoodI (l2000ll . Since this library includes 
only cool stars, they added the oretical hot star s. A simi¬ 
lar approach was carried out by iMarastonI (l2005l l who also 
took advantage of the empirical spectra of the library by 
iLancon fc_Wood (200(]|) f or coo l stars. The set of models de- 
rived bv iMaraston et al] (l2009l ~) is comp letely based o n the 
empirical low resolution stellar library bv |Pickl3 (Il998fl . Re¬ 
cently, Meneses-Goytia et al. (2014) took advantage of the 
IRTF-library to create SSP models ranging up to 2.4/rm. 
Their approach is quite similar to the one conducted by us, 
since the input stellar library and also the routines used to 
create the model spectra are identical. Here, however, we 
focus on different wavelength ranges and pursue different 
global aims (details see below). 

The main goal of our work as presented in this paper 
is to extend the SSP models developed in the optical 
spectral range between A = 3465 A and A = 9469 A by 
IVazdekis et al.l (l2003l . |2010| . I 2 OI 2 I I to the IR wavelength 
regime. The IR spectral range offers certain advantages 
compared to the optical. Firstly, it is much less affected by 
dust extinction. Moreover, IR light traces cool, old stars 
which dominate the baryonic mass in galaxies. Finally, 
these wavelengths are also very suitable to quantify the 
contribution of asymptotic giant branch (AGB) stars. 

Our models are based on the same modelling routines (see 
Section [8l) as the ones d escribed in the aforementioned 
papers by IVazdekis et aP (I 200 II 2 OI 0 I. I 2 OI 2 II . As described 


there, we followed an approach which is as empirical as 
possible. Hence, contrary to other authors, our SSP models 
are - apart from two gap regions, see Section!?]- completely 
based on an empirical stellar library. In addition, we also 
made use of extensive photometric libraries rather than 
theoretical stellar atmospheres in order to transform the 
theoretical parameters of the isochrones to the observational 
plane. Furthermore, we did not enlarge our empirical stellar 
library by including artificial stars. Our models focus on the 
spectral range between 2.5 /rm and 5 /rm in which the Spitzer 
[3.6] —[4.5] color can be found. However, in a forthcoming 
paper we will combine our m odels with the already existing 
ones of (IVazdekis et al.ir2012l l. This will enable us to provide 
the first empirical homogeneous stellar population synthesis 
models ranging from the optical to the IR. 

The paper is structured as follows. In Section 2, we present 
the IRTF library upon which our SSP models are based. 
After that, we describe all the steps which had to be car¬ 
ried out in order to prepare the IRTF stars for its use in 
SSP modelling. We had to test whether the library stars 
had been correctly flux-calibrated (see Section 3), to find 
out about the spectral resolution of the library (Section 4) 
and to determine as accurately as possible the stellar at¬ 
mospheric parameters of all of the stars (see Section 5). In 
addition, we checked all of the library stars for peculiarities 
(see Section 6) before finally correcting the chosen stars by 
fitting gaps occurring due to telluric absorption lines and ex¬ 
trapolating them to a common dispersion over all the wave¬ 
length range (see Section 7). In Section 8 we explain the 
modelling process, the main ingredients used and in which 
parameter range our models can be safely applied. We study 
the behaviour of colors measured from our model spectra 
as a function of age and metallicity and compare them to 
photometric predictions from MILES llVazdekis et al.l[2010h 
(see Section 9). In Section 10, we compare our SSP mod¬ 
els to other models available from the literature. Section 11 
is dedicated to the discussion of mass-to-light [M/L) ratios 
measured from our models and to the investigation of the 
dependence of our models on the IMF. Finally, we compare 
the colors deduced from our SSP models to observed colours 
of globular clusters and of early-type galaxies (Section 12) 
and end the paper with a discussion and some concluding 
remarks (Section 13). 


2 THE IRTF LIBRARY 

The IRTF spec tral library (ICushing et al.l l2005l : 
iRavner et al.l l2009l 'l contains the spectra of 210 cool 
stars in the NIR and IR wavelength ranges. All of the stars 
were at least observed in the spectral range from 0.8 - 
2.4 ^m, but a significant fraction extends up to 4 pm, and 
about half of the stars even up to 5 pm. 

The stars were obser ved with the med ium-resolution 
IR spectrograph SpeX llRavner et al.l [2003li at the 3.0 m 
IRTF on Manna Kea, Hawaii, with a resolving power of 
R = A/AA «2000. The total wavelength range up to 
5 pm was covered in only two cross-dispersed instrumental 
settings. The signal-to-noise ratio (S/N) is better than 
« 100 except for regions of poor atmospheric transmission 
and for A > 4 pm and all stars have been observed at an 
airmass < 2 for good telluric corrections. 
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The observed sample is composed of stars of spectral types 
F, G, K, M, some AGB-, carbon- and S-stars, of luminos¬ 
ity classes I-V. Since it was aimed to achieve a high S/N- 
ratio also in the IR (2.4 - 5pm), the selected stars are all 
relatively bright and situated in the Solar neighbourhood. 
Consequently, most of the stars of the library are of So¬ 
lar composition. Unfortunately, this caveat will prevent us 
from constructing SSP models of highly sub-and supersolar 
metallicities (see Section [8l. 


3 TESTING THE FLUX CALIBRATION OF 
THE LIBRARY STARS 

Since an important application of our SSP models will be 
to extract colours from the calculated model spectra and 
compare them to observations, we have to assess the flux 
calibration of the stellar spectra. Thus, when integrating 
the stellar spectra, each star needs to have the correct flux 
at a given wavelength - at least with respect to the fluxes at 
other wavelengths. Modelling requires a correct relative flux 
calibration, but absolute fl ux calibration i s not a necessary 
prerequisite. According to iRavner et al.l (l2009l l. the IRTF 
spectra were relatively flux-calibrated and at the same time 
corrected for absorption due to the Earth’s atmosphere by 
dividing them by a telluric correction spectrum constructed 
from a theoretical model spectrum of Vega scaled to the 
spectrum of an observed nearby standard star. The absolute 
flux-calibration was done by using the Two Micron All Sky 
Survey f2MASS. ISkrutskie et al.]l2n0flh JHKs photometry. 
To test the flux calibration qua lity of the IRTF spectra car¬ 
ried out bv iRavner et al.l (120091 '). we calculated the 2MASS- 
colors (J — H), (J — Ks) and {H — Ks) from the spectra 
by integrating over the spectr al flux in the 2MASS fil ters 
and using a Vega spectrum of ICastelli fc Kurucd lll994l 'l as 
reference. Then we plotted our calculated colours versus the 
respective colours extracted from the 2MASS catalogue. In 
general, we observe quite a big scatter around the identity 
line, particularly for {H — K). 

As a next step, we created a colour-colour plot of (J — K) 
versus (J — H) (see Figure [T]). In this plot, we did not only 
display the respective colours for the IRTF stars as calcu¬ 
lated from their spectra and extracted from the 2MASS cat¬ 
alogue, bu t we als o inclu ded the colours of the stars of the 
library by IPicklesI l|l998h , once obtained by us from their 
spectra and again catalogue d by the latter, an d additionally 
the JHK-colouis derived by iKoornnee 3 (Il983ll . Figure [U re¬ 
veals that the colours of all these samples agree well amongst 
each other, apart from a significant number of colour val¬ 
ues for the IRTF stars originating from the 2MASS cata¬ 
logue. Consequently, we conclude that in many cases the 
colours for the IRTF stars deduced from the 2MASS cata¬ 
logue are somewhat questionable since they seem to be in¬ 
accurately dete rmined. This picture of 2MASS photometry 
is confirmed by iTurneil (120101 . [2oT3) who claim that it is of 
uneven quality and low intrinsic precision with errors larger 
than the cited values for most of the stars in the Galactic 
plane. However, the fact that the colours which we deter¬ 
mined directly by integrating over the IRTF spectra are in 
good agreem ent with the corresponding co lours from other 
sources fe.g.. IKoornnee 3 l 19831 : [Pickleslfl998l l proves that the 
flux calibration has been done correctly for the stars of the 



Figure 1. (J — K) versus (J — H) colour-colour plot. Black crosses 
indicate colours which we calculated from the IRTF spectra, blue 
stars stand for colours of the IRTF stars extracted from the 
2MASS catal ogue, whereas red diamonds delineate the colours 
of the stars of lPickl^ lll998l ) as catalogue d there. Colou rs which 
we obtained directly from the spectra of IPickl^ lll998ll can be 
identified by green triangles a nd magenta squa res illustrate the 
corresponding JHK-colours of iKoornneell (ll983l V 


IRTF library. This is obviously the most important result of 
our testing. 


4 DETERMINING THE SPECTRAL 
RESOLUTION 

It is important to find out the spectral resolution of the 
IRTF stars since this will be also the one of our SSP models. 
Furthermore, we need to know how the spectra have to be 
treated during the modelling process. We have to determine 
the resolution as a function of the wavelength A. Equation 
[T] shows the relation between the respective parameters: 


FWHM a • 2.35 


A 


1 

R 


( 1 ) 


In Equation (ID, R deno t es the resolving power which, ac¬ 
cording to iRwner et al.l ll2009ll . is « 2000. We determined 
the spectral resolution by fitting a combination of spectra 
of the most recent relea se of the theoretical Phoenix li¬ 
brary (lAllard et al.ll20f^ l to the IRTF spectra. This library 
comprises spectra of a large parameter range, but we only 
used the model spectra extending from log(g) = —0.5 to 
log{g) = 5 and Teg = 5200 K to 6000 K since we carried 
out this fitting only for IRTF stars within this temperature 
range. We finally recovered the FWHM and the radial ve¬ 
locity of our set of tes t stars by applying the penaliz ed pixel 
fitting method ('pp xf. ICappellari fc Emsellem 2004l l follow¬ 
ing the approach of iFalcon-Barroso et al.l ( 201lf) . The spec¬ 
tral fitting was done separately for seven spectral ranges of 
always around 2000 A width covering the range of the test 
spectra until the end of the K band. Since the theoretical 
spectra of the Phoenix library do not model very accurately 
the spectral features beyond the K band, we just give a 
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Figure 2. Upper panel: Dependence of the FWHM on wave¬ 
length. Shown are the average values for the respective wavelength 
bins obtained from the sample of 18 test stars of temperatures be¬ 
tween 5200 and 6000 K. The red errorbars indicate the errors of 
the mean values. The black line indicates the FWHM correspond¬ 
ing to a resolving power of R = 2000 over this wavelength range. 
Lower panel: Dependence of the velocity dispersion cr on wave¬ 
length. The black line delineates the best fit to the data points. 


single resolution estimate for the reddest spectral range. Fi¬ 
nally, we calculated one average value for all of our eight 
wavelength bins from the 18 test spectra which we selected. 
The upper panel of Figure[5]shows that the FWHM increases 
with wavelength from around 5 A in the J band to around 
10 A in the K band. Moreover, it looks like this trend con¬ 
tinues beyond the K band. However, we should be careful 
drawing such a conclusion given the fact that we were only 
able to obtain one further value for the whole wavelength 
range between the end of the K band and 5 /rm. 

As expected from Equation [T] the lower panel of Figure 
[2] shows that a is virtually constant throughout the whole 
spectral range at a value of ~ 60 km/s. A fit to our data fur¬ 
nishes the exact value of ct = (60.6 ± 1.1) kms“^ (see Figure 
Such a value is appropriate to study most galaxies. 


5 DETERMINATION OF THE STELLAR 
ATMOSPHERIC PARAMETERS 


To use the IRTF stars in SSP modelling it is crucial to de¬ 
termine their stellar atmospheric parameters effective tem¬ 
perature Teff, surface gravity log(g) and metallicity [Fe/H] 
as accurately as possible. To measure these parameters, we 
applied two different methods. The first one consists of 
compiling the values from the literature and subsequently 
adjusting them to a master template follo wing the ap¬ 
proach carried out by ICenarro et al.l (120071 ') (see Subsec¬ 
tion I5.1|l . As a second method, we adopted a number o f 
colour-te mperature relations from lAlons o et al.l (Il996l . [l999l l 
and from IConrov fc van Dokkumn 20121 ') as well as the re¬ 


lation between temperat ure, bolometric magn i tude, mass 
and surface gravity from IStraizvs Sz Kuriliend (Il 98ll) fsee 
Subsection 15.21) . The relations bv lAlonso et alTl|l996l . ll999l ') 
are the same ones which we later employ in our modelling 
procedure to transform the theoretical parameters of the 


isochrones to fluxes and colours (see Section l8.ll) . Finally, 
we checked whether our stellar atmospheric parameters ob¬ 
tained with these two methods were in accordance with 
the ones expected for a star of the respective spect ral type 
and lu minosity class according to the catalogue by IPicklea 
l|l998h . Moreover, we also cross-checked our final set of de- 
rived atmosp h eric p arameters with the sample compiled by 
ICesetti et al.l ll2013l l from the literature for many of the 
IRTF stars. 


5.1 Compiling the stellar atmospheric parameters 
from literature 

We started by querying the literature in order to find all rel¬ 
evant sources containing spectroscopically determined stel¬ 
lar atmospheric parameters. Then, we checked whether any 
of the 180 stars that we used from the IRTF library were 
contained in the different references. Whenever this was the 
case, we matched the respective sample of stars from the 
literature with the 985 stars of the MILES stellar library 
dSanchez-Blazauez et al.l 120061 ') to find stars in common. In 
order to be fully consistent with our model predictions in 
the optical range, the stellar parameter s of the MILES stars 
as determined bv ICenarro et aO (l2007l l served as the basic 
reference sample against which all the other references were 
calibrated. Finally, we ended up with 17 sources from the lit¬ 
erature which have stars in common with both the MILES 
and the IRTF stellar libraries (see Table [T|. 

After that, the three parameters Tefj, log( 5 ) and [Fe/H] 
of the stars in common between one literature source and 
MILES were separately compared to each other. For all 
three parameters, we subsequently obtained the best-fitting 
lines mapping the relations between the MILES values and 
their equivalents from the respective literature reference. 
When carrying out these fits, we first used all stars in 
common between MILES and the literature source. In the 
next step, we discarded all those stars showing a value 
A(parameter) = jparameterMiLES — parameter reference j 
larger than la (i.e., 1-a clipping) and performed another 
fit. The offset parameter A and regression coefficient B of 
this latter fit were then used to correct the stellar parame¬ 
ters of all the stars of the respective reference. Like this, we 
also enlarged significantly the initial reference sample which 
consisted of MILES stars only. Consequently, also those 
references (Santos et al.|[2004|, |20 05l:ISousa et al.||2008l, |201l|; 


KxMvul^ 20071 : Kovtvukh et al.l feoosi : Lvubimkov et al.l 

2 OIOII which a priori did not exhibit enough stars in common 
with the MILES library to obtain reliable fits could get 
adjust ed to the MILES scale. According to ICenarro et al.l 
(( 20 ^, for this a minimum number of 25 stars in common 
is needed. 

For all of our references, we also created plots showing the 
differences between the respective parameter values and the 
MILES ones versus the MILES values of all three different 
parameters. In the case of most of our references, the a stan¬ 
dard deviations of the parameter differences do not exceed 
75K in Teff, 0.2 in log(( 7 ) and 0.05 in [Fe/H] (see Figure jS]). 
Finally, we found that the parameters of 44 of our 180 IRTF 
stars could be directly adopted from MILES. 32 further stars 
were contained in at least one reference, whereas 20 addi¬ 
tional ones were discovered in several references. The final 
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Table 1. Literature sources containing stellar atmospheric parameters of stars in common with the MILES and IRTF libraries. 


Source 


IRTF library 
^Prugnie^^t_^ ^2007) 

Wi^^t_^ (^20]d ) 

_Cenarro_e^_^ (_200^ [T^g, log(g)] 
_Cenarro_et_^ ^200^ [Fe/H] 
_Gr^_e^_^ (^200^) 

Va]enti_&_^ischer ^200^ 
IIekker_&_^elendez (200^ 

Takeda (20071 

^Mishenina_e^_^ (2004 , ^006 , ^00^ 
I^ovt^uM^_e^_^ (2003 , ^00^ ["^eff] 
I<ovt^ukl^_e^_^ (2008) [log(g), [Fe/H]] 
^Lucl^_&_JIeiter (2007) 

_Santos_e^^l^ (2004, ^00^ 

_Sousa_et_^l^ (2008, ^01^ 

I^ovt^uMi (2007) 

^L^ubimkOT^t_^ (201^ 
^Edrardssor^_et_a^ (1993 ) 
^Borkova_&_^arsakov (2005) 

Thevenin (19981 


stars in common 
with the MILES library 

47 

302 

373 

348 

308 

43 

64 

72 

90 

149 

92 

27 

53 

35 

37 

11 

5 

94 

210 

273 


stars in common 
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Figure 3. Differences in temperature between MILES and two 
literature sources as a function of temperature (left-hand pan¬ 
els) and metallicity (right-h and panels), resp e ctively. These two 
sources are the catalogue of iMishe nina et ^ ll2004l . |2006|. [2008^ 
(upper panels) and the one of iTheveni nl ll 19981) (lower panels). 
The red solid lines are the ones indicating the best fits, the black 
dashed lines represent the mean temperature differences between 
the two samples and MILES, whereas the black dotted lines show 
the one cr standard deviations from these mean differences in tem¬ 
perature. 


parameters p of the latter on es were calcula t ed ap plying the 
following equation, given bv ICenarro et all (|2007l) : 

N 

E Pi/^i 

p — - .S APscolortemperaturerelations (2) 

E 

i=l 

Here, p* denotes the parameter of one literature reference 


corrected to the MILES scale, and at corresponds to the 
standard deviation of this reference with respect to the 
MILES system. 

We were not able to derive stellar atmospheric parameters 
for 79 of our stars from the literature. For these stars, colour- 
temperature relations had to be applied (see next Section 

EU. 


5.2 Applying colour-temperature relations 


Since literature values are not available for all of our used 
IRTF stars, we also followed a second, different approach to 
determine the effective temperatures of the stars of the IRTF 
library. Theref ore, we used the colou r-temperature relations 
deter mined bv lA lonso e t al.l (Il996l) for dwarfs and for gi¬ 
ants llAlonso et al.l 19991 ) of spect ral types F0-K5 for colours 
ranging from optical to the NIR. lAlonso et al.l lll996l . Il999l) 
indicated in which colour ranges their relations are appli¬ 
cable. These are metallicity-dependent relations, although 
for the (J — K) colour the sensitivity to this parameter is 
negligible. We consequently calculated the J and K mag¬ 
nitudes from the spectra of all of our stars by convolving 
them with the respective filter response functions. For the 
stars of the IRTF library which only extend up to 2.4 pm 
we had to use the slightly metal-dependent relation for the 
(J — H) colour in order to calculate the temperatures, since 
their spectra do not cover the whole spectral range of the 
K filter . It should b e noted that the J, H and K filters 
used by lAlonso et ahl (Il996l . Il999l ) for their observations at 
the Carlos Sanchez Telescope (TCS) at the Observatorio del 
Teide in Tenerife are non-standard on es. The filter r e spons e 
functions are give n in the appendix oflAlonso et al.l ll 19941 ). 
The relations of lAlonso et al.l lll99(j . 19991) are however 
not valid for M-stars. For them, we adopted the t e mper - 
ature values calculated by IConrov fc van Dokkurnl (|2012|) 
for these stars using the temperature-colour-relations by 
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Figure 4. Differences b et ween e ffectiv e t emperatures determined 
following lAlonso et alj lll996l . Il999tl ; IConrov fc_vBJi_^>Qkku^ 
j2^^) and according to a literature compilation Jcenarr^^^^ 
l2007h versus the latter ones. The black pluses indicate Teff val- 
ues of the respective stars directly from the MILES catalogue 
(|Cenarro et al.ll20Q7h . whereas the blue asterics and green dia¬ 
monds indicate that the corresponding Teff is based on one or 
several literature sources, respectively. The solid black line marks 
the identity between the e ffect ive temperatures dete rmined using 
lAlonso et al.l lll996l , Il999ll and ICenarro et al.l ll2007t) . 



Figure 5. Distribution of our IRTF stars in the log(^)-Teff 
parameter space. Blue pluses indicate stars of a metallicity 
lower than [Fe/H] = —0.6, orange stars stand for stars in the 
metallicity range between [Fe/H] = —0.6dex and [Fe/H] = —0.3, 
whereas red diamonds delineate stars of metallicities between 
[Fe/H] = —0.3 and [Fe/H] = —0.1. Stars of around solar metallic¬ 
ities (—0.1 < [Fe/H] < 0.1) can be identified by green triangles 
and black squares illustrate all stars of metallicities greater than 
[Fe/H] = 0.1. 


iRidgway et al.l (fl980l) andjPerrin et al.l (|l998l) for M giants 
and bv fcasagrande et aP (120081 . l2010l ) for M dwarfs. 

Figure U compares the temp eratures calc u lated using 
the aforementioned relations of Alonso et al.l (|l996l . 119991 1 
and I Conroy &: van DokkunJ lj2012h to the temperatures 
determined b e foreh and by means of the approach of 
ICenarro et alJ ll2007h for those stars with atmospheric pa¬ 
rameters in the literature. 

The variations are reasonable in general, although they are 
quite pronounced for some individual stars and the temper¬ 
atures from the color-temperature relations appear to be 
a bit lower on average than the ones from the literature 
compilation. Consequently, in order to decide which of the 
two temperature valu es to adopt we performed a cross-check 
with the catalogue bv IPickl^ (Il998l b This library consists of 
131 stellar spectra extending from 1150 A to 10620 A (half 
of them even further until 25000 A) which were artificially 
compiled based on several sources. We searched for coun¬ 
terparts having the same spectral types like our IRTF stars 
among the Pickles stars. Before, we had a detailed look at 
all the spectra by eye in order to make sure that their spec¬ 
tral types were properly determined. Then we checked what 
temperature a star o f a sp ecific spectral type should have 
according to |Pickl3 (Il998fl . in order to get an idea about 
which of our two temperature values is most likely the more 
reliable one. In the end, we substi tuted the temperature val¬ 
ues from our ICenarro et ahl ll2007l '> approach by the ones cal- 
c ulated using the colour-tem perature(-metallicity) relations 
of lAlonso et al.l lll996l . fl999ll for a number of stars. For most 
of these stars, only one literature reference, different from 
MILES, had been found. For a small number of stars, we 
used the fast approach of taking directly the temperature 
values given in the Pickles-library for a star of a similar 
spectral type. 


For the five carbon stars in our sampl e, we adopted the ef - 
fective temperature values obtained bv iBergeat et al.] (1200 ill 
as the arithmetic mean of the temperatures estimated from 
the integrated flux and from calibrated colour indices. 

We also calculated log( 5 r) for all IRTF stars for which 
this was feasible due to their spectral type and luminos¬ 
ity class by applying the relation between surface grav- 
ity, temperature, stel l ar ma ss and bolometric magnitude of 
IStraizvs fc Kuriliend ([TosJ). For those stars for which this 
calculation yielded a value of the surface gravity which is 
in better agreement with the value expected from the lumi¬ 
nosity class given in the IRTF library, we substituted the 
value. 

Table in the appendix contains our final set of stellar 
atmospheric parameters as used later in the SSP modelling 
(see Section [8]). In this Table, we also indicate how the pa- 
rameters were determ ined. In general, the values from our 
ICenarro et al.l ll200lt) approach were always preferred over 
the ones from the colour-temperature relations as long as 
they wer e in reasonab le agreement with the expected val¬ 
ues from IPicklesI lll998ll . Figure [S] illustrates the distribution 
of these parameters in the log(g)-Teff plane. The respective 
stars are colour-coded according to their metallicities. From 
Figure [S] it is obvious that a good parameter coverage in 
Tefi and log( 5 r) only exists for stars of around Solar and for 
slightly sub- and supersolar metallicities. Subsequently, we 
model SSPs just within this metallicity range (see Section 

B- 

6 STELLAR CONTENT OF OUR MODELS 

In this Section, we discuss some aspects of the final sample 
of stars on which our SSP models are based. At first, we 
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comment on the stars which we hnally had to discard or at 
least reduce in weight due to various reasons (see Subsection 
16.111 . Since the SSP spectra in the NIR and IR depend mostly 
on the coolest stars, in Subsection we describe what 
kind of cool stars entered the models. Moreover, we give 
the fractions of the final modelled bolometric flux due to a 
particular subsample of cool stars. 


6.1 Checking for peculiar stars 

For their later use in SSP modelling, it is very important 
to make sure that we only consider spectra of normal stars 
- i.e., stars without any peculiar spectral features such as 
very strong emission lines. In order to decide which stars to 
discard from our sample of IRTF stars or of which we should 
at least reduce the weight with which they will contribute 
to the modelling process, we took advantage of the so-called 
interpolator (se e Section J82L f or a d etailed explanation see 
the appendix of IVazdekis et al]l2003l l. 

We then used this programme in order to create artihcial 
stars with the same parameters as all the stars of the IRTF 
library. When carrying out this interpolation, we always ex¬ 
cluded the actual star whose spectrum should be reproduced 
artificially from the list of stars contributing to the inter¬ 
polation. Then, we created plots showing the residuals be¬ 
tween one spectrum and its artificial counterpart. Before 
analysing these residual plots statistically, we masked out 
some wavelength regions hampered by strong telluric fea¬ 
tures and spikes (see Figure li). Stars showing a standard 
deviation smaller than around 0.03 and a maximum and 
minimum value deviating by no more than around 20 per 
cent from unity qualihed to contribute with their full weight 
to our models. We had a detailed look at all those stars 
which exceeded these limits. For many of these stars, the in¬ 
terpolator was simply unable to create a reliable counterpart 
of the same parameters due to the limited number of stars 
in the respective regions of the parameter space. This was 
in particularly the case for stars of extremely high and low 
metallicity, C-stars, S-stars, dwarfs of very low temperature, 
dwarves of high temperature and comparably low surface 
gravity, high-temperature giant stars and low-temperature 
supergiants. Whenever it was feasible, we rather reduced the 
weights with which these stars contribute to the SSP mod¬ 
els than discarded them completely. The exceptions were 
star HD 142143) which features a very bad residual plot 
with a never constant baseline, stars HD 14386 and HD 
69243 because of their multiple emission lines, and stars GJ 
1111 and HD 6474 for which wrong temperature values were 
clearly assumed. Moreover, the four S-stars (BD-l-442267, 
HD 64332, HD 44544, HD 62164) had to be eliminated, as 
they are long-period variables and also contain some pecu¬ 
liar spectroscopic features. We reduced the weight of 18 stars 
to 50 per cent of the weight given to a good star (HD 40535, 
HD 87822, HD 17918, HD 126660, HD 202314, HD 179870, 
HD 212466, BD+60265, HD 339034, HD 39801, HD 206936, 
HD 95735, HD 14488, G1 268, HD 196610, HD 42581, HD 
35661, HD 14404) and of one star (G1 406) to only 25 per 
cent of the weight of a good star, and excluded a further 
nine stars completely (see above). 



Figure 6. Example of a residual plot. We divided the spectrum 
of HD 21018 (Teff = 5353 K, log(g) = 3.07, [Fe/H] = 0) by the 
one of its artificially created counterpart obtained by using the 
interpolator. The red line runs parallel to the x-axis at j/ = 1. 
The grey-shaded areas are the parts of the spectrum which were 
masked out when getting the statistics (mean value, standard de¬ 
viation, maximum and minimum value) of the plot, since they are 
more affected by noise and telluric features and thus less reliable 
for these purposes (see text). 


6.2 Cool stars 

Figure [5] in Section reveals that despite the generally 
poor coverage of the stellar atmospheric parameter space, 
potentially very important and relevant stars for SSP mod¬ 
elling in the IR like cool dwarfs and supergiants form an 
integral part of our finally used stellar sample. Hence, our 
sample of 180 stars contains in total 32 AGB stars, amongst 
them 16 very cool supergiants of temperatures between 
Teff = 2300 K and Tefi = 3500 K and 17 cool dwarfs with 
temperatures ranging from Teff = 2700 K to Tgff = 3840 K. 
Table [2] lists all the mentioned cool dwarfs and supergiants 
together with their stellar atmospheric parameters. 
Consequently, we are able to reach masses as small as 
O.15M0. In the standard case of a Kroupa-like IMF, the con¬ 
tribution of these cool very low main sequence red dwarfs of 
masses smaller than 0.5Mq to the total bolometric flux of 
one particular SSP-model is between 1 % at ages of around 
1-2 Gyr and 5 — 6 % at the oldest ages, irrespective of 
the metallicity. For a more bottom-heavy IMF (see Section 
EJI, these fractions increase to values of between 6 and 
14 %. As expected, the cool AGB stars are particularly im- 
por tant at the lower ages between 1 and 2.5 Gyr (compare 
also lMaraston|[2005l : lMaraston et al.ll2009ll . In this age range, 
the AGB stars are responsible for more than 20 % of the 
total bolometric flux of the models, whereas this percent¬ 
age has dropped down to only « 10 % at the oldest ages 
of 10 — 12 Gyr. Among the AGB stars of our sample, also 
some carbon stars can be found. Our input sample of 180 
stars contains five carbon stars, stars HD 31996, HD 57160, 
HD 76221, HD 44984 and HD 92055 of temperatures be¬ 
tween 2300 and 3300 K. Carbon stars contribute up to 5 % 
of the total bolometric flux in the case of solar and higher 
metallicities at the youngest ages covered here (1 — 2Gyr). 
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This fraction increases for subsolar metallicities to values of 
« 10 %. For older ages, carbon stars do not play a role and 
their contribution to the total bolometric flux is negligible. 
Oxygen-rich stars, however, are not included in the IRTF 
library or at least it is is not explicitly pointed out. All the 
presented estimates of the percentages of the total bolomet¬ 
ric flux assigned to the various ca tegories of cool stars are 
based on the BaSTI-isochrones of IPietrinferni et al.l (l2004l f 
(more details on the modelling process and on the necessary 
ingred ients in Section [51 ) . Alth ough it is not stated explic¬ 
itly in IPietrinferni et al.l (l2004l . they do take the so-called 
third dredge-up, during which convection transports the car¬ 
bon produced in the stellar interior by helium fusion up to 
the stellar surface, into account ( Cassisi, private communi¬ 
cation and IVazdekis et al.l (120151 11. For doing so, they ap¬ 
ply the analytical relationships for the trend of the surface 
luminosity and effective t emperature as a function of th e 
CO core mass provided bvlWagenhuber fc TuchmanI lll996ll : 
IWagenhuber fc GroenewegenI (119980 and a correction to the 
(J — K) colour whi c h dep ends on the effective temperature 
from iBergeat et ^ (1200 ih . 


7 PREPARING THE STELLAR SPECTRA 


In this Section, we will briefly describe the final treatments 
which we had to carry out on the spectra before being able 
to use them for SSP modelling. First of all, it turned out 
that all of the IRTF spectra show three characteristic gaps, 
one between 1.80 and 1.87 /rm (i.e., between the J and the 
H band), another one between « 2.50 and 2.80/rm and a 
third one between « 4.20 and 4.50 fim. The exact limits of 
the latter tw o ones vary slightly f rom spectrum to spectrum. 
According to lRavner et al.l (l2009tl . the first gap is mainly due 
to the instrumental setting, whereas the second and third 
ones are caused by strong telluric absorption in the respec¬ 
tive wavelength ranges. We interpolated the spectra linearly 
within the first two gaps using IRAF. Whenever it was feasi¬ 
ble, we replaced th ose two gaps by the theoretical spectra of 
the BaSeL-library jLeieune et al.lToOTl . 1 199^ : 1 Westera et al.l 
I 2 OO 2 I I enabling us to include also some spectral features. 
The spectra of this library extend over the large wavelength 
range 91 A to 1600000 A, however at a rather low resolution 
of only 10 - 20 A, and were compiled for a close grid of stellar 
atmospheric parameters. The upper two panels of Figure [7] 
show the spectral ranges around the first two gaps of the star 
HD 213306 (Tefi = 6051 K, log(g) = 1.69, [Fe/H] = -0.12) 
containing also the corrections carried out within these gaps, 
based on the corresponding BaSeL-spectrum. 

There is a strong molecular CO absorption band in the range 
betw een 4.20 /rm and 5 /xm which affects the 4 .50 /rm filter 
(e.g., iMonson et al.l I2OI2I : IPeletier et al.l I2OI2I : iMeidt et al.l 
l2014h . Unfortunately, this CO absorption falls in the gap 
present in all of the IRTF spectra in the region between 
4.2 jj,m and 4.5 /rm, and thus prevents us from properly ac¬ 
counting for this flux loss in the 4.50 /rm band. As can be 
seen from Figure 10 in iMonson et al.l (l2012fl . this CO absorp¬ 
tion affects stars of temperatures cooler than about 4500 K, 
i.e., the ones which contribut e most to our models in the IR 
(see also IPeletier et ^l2012l l. 

In order to correct our spectra for the effect of the CO in 
the 4.5 pm band, we made use of the theoretical Phoenix 


Table 2. List of the coolest stars of the library: The first column 
gives the effective temperatures Columns 2 and 3 contain the 
surface gravities log(g) and the metallicities [Fe/H], respectively. 
The official star name can be found in Column 4 and the type of 
star (carbon, AGB or M dwarf) is indicated in Column 5. 


Teft [K] 

log(5) [dex] 

[Fe/H] [dex] 

star name 

type 

2314 

-0.56 

0.260 

HD31996 

carbon 

2645 

-0.45 

-0.039 

HD76221 

carbon 

2698 

4.63 

-0.149 

G1406 

M dwarf 

2750 

-0.56 

-2.609 

HD207076 

AGB 

2794 

5.35 

-0.147 

G151 

M dwarf 

2951 

5.27 

-0.151 

G1866 

M dwarf 

2965 

-0.56 

-0.060 

HD92055 

carbon 

3030 

5.20 

-0.544 

G1213 

M dwarf 

3035 

-0.45 

-0.074 

HD44984 

carbon 

3060 

0.31 

-1.016 

HD196610 

AGB 

3108 

5.00 

-0.247 

G1273 

M dwarf 

3111 

5.22 

-0.718 

G1299 

M dwarf 

3197 

4.57 

0.000 

G1581 

M dwarf 

3236 

0.73 

-2.609 

HD156014 

AGB 

3250 

0.30 

-1.005 

HD18191 

AGB 

3278 

-0.25 

-0.069 

HD339034 

AGB 

3295 

-0.45 

0.046 

HD57160 

carbon 

3304 

4.87 

-0.158 

G1388 

M dwarf 

3319 

4.63 

0.000 

G1381 

M dwarf 

3349 

0.47 

0.026 

HD40239 

AGB 

3397 

0.90 

-1.396 

HD214665 

AGB 

3420 

0.50 

0.140 

HD175865 

AGB 

3428 

1.00 

-1.464 

HD204585 

AGB 

3434 

0.12 

-2.267 

HD108849 

AGB 

3454 

0.37 

0.183 

BD-l-60265 

AGB 

3509 

-0.16 

-0.772 

HD 14488 

AGB 

3544 

0.98 

-0.959 

HD19058 

AGB 

3550 

0.04 

0.030 

HD39801 

AGB 

3551 

-0.16 

-0.693 

HD14469 

AGB 

3571 

0.99 

-0.508 

HD4408 

AGB 

3576 

0.76 

-2.609 

HD94705 

AGB 

3594 

0.52 

0.177 

HD236697 

AGB 

3611 

4.71 

-0.192 

G1806 

M dwarf 

3617 

0.06 

-0.201 

HD35601 

AGB 

3640 

0.32 

-0.070 

HD 14404 

AGB 

3650 

4.69 

-0.149 

HD209290 

M dwarf 

3730 

0.90 

0.270 

HD219734 

AGB 

3731 

4.65 

-0.161 

HD42581 

M dwarf 

3737 

4.90 

-1.500 

HD36395 

M dwarf 

3746 

-0.04 

0.046 

HD206936 

AGB 

3750 

4.75 

-0.156 

G1268 

M dwarf 

3781 

0.50 

-0.458 

HD 10465 

AGB 

3782 

0.74 

-0.928 

HD216946 

AGB 

3797 

0.00 

0.064 

HD212466 

AGB 

3828 

4.90 

-0.200 

HD95735 

M dwarf 

3843 

4.49 

-0.408 

HD201092 

M dwarf 

3848 

0.27 

-0.151 

HD185622 

AGB 

3863 

0.78 

-0.626 

HD23475 

AGB 


library (lAllard et al.ll20l^') for filling in the gap. Thanks to 
their high resolution of 0.1 A in our IR wavelength range, 
the spectra of this library account for the effect of the CO 
absorption. In addition, they cover all of our temperature 
range from « 2000 K to « 7500 K with a spacing of 100 K, 
contain surface gravities ranging from —0.5 dex to 5.50 dex 
in steps of 0.5 dex and are available for metallicities between 
[Fe/H] = —4.0 dex and [Fe/H] = 0.5 dex with A = 0.5 dex. 
The lower panel of Figure [7] illustrates the fitting of this 
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X [|im] 

Figure 7. The three ranges of the spectrum of HD 213306 show¬ 
ing gaps. The original spectrum is delineated with a black dotted 
line. We overplotted the gap-corrections in red. 


8 INFRARED SSP MODELS 


SSP models are a very important tool for studying clus¬ 
ters and galaxies and in particular to determine their ages 
and metallicities. These comparisons ca n be carried out by 
means of full spectrum fitting (see, e.g.. lKoleva et ani201l[l 
or by using information about integrated colours or line 
strength indices. 


SSP models are readily available (e.g., 

Vazdekid 19991: 

Bruzual & Charlotll2003l: iGonrov et al.ll201C 

: Vazdekis et al. 

201 c 

: Maraston & StrombackI 201ll: Conrov & van Dokkum 

2012 

: Vazdekis et al. 2012|). whereas in the NIR and IR 


wavelength ranges, their number is more limited. This 


is particularly true for the IR. Whereas examples fo r 
the NIR includ e for e xample iMouhcine fc LanconI ll2002ll : 
iMaraston et al.l ll2009ll . there do not exist any SSP models 
beyond the K band which are mainly based on empirical 
stellar spectra and which predict model spectra at moder¬ 
ately high resolution. 


third gap by the theoretical Phoenix spectrum with the clos¬ 
est agreement in the stellar atmospheric parameters to the 
ones of the respective IRTF star. Finally, this approach of 
correcting for the CO absorption feature enables the reader 
to deduce colours involving the 4.5 /rm band from our SSP 
model spectra without having to apply a missing-flux cor¬ 
rection accounting for the gap. 

We note that the spectra of the IRTF library are not given 
with a common pixel size, but that it grows from around 
2Apixel“^ at 0.8/rm to around 10Apixel“^ at 5/rm. In or¬ 
der to use the spectra for SSP modelling, it is necessary to 
have them at a common pixel size. We thus rebinned them 
to AA = 2.73 A, since this is the value corresponding to a 
wavelength of around 1 /rm. 

We still had to extrapolate the spectra of all of our used 
IRTF stars which do not already intrinsically extend up to 
5 fj,m to this value. In order to fulfill this task, we made again 
use of the interpolator. At hrst, we created artificial coun¬ 
terparts reaching up to 5 /rm for all of the IRTF stars whose 
spectra did not extend beyond 4 pm. This interpolation was 
based on all stars with spectra observed up to 5 pm. Then, 
we combined the artificially created part of the spectrum 
between 4 and 5 pm with the actual spectrum of the respec¬ 
tive star reaching up to 4 pm. Therefore, we selected the last 
part of the original stellar spectrum from around 39500 to 
40000 A and calculated the mean offset in this spectral range 
between the spectrum of the 4 pm IRTF-star and the corre¬ 
sponding star of the same parameters and added this offset 
to the latter one. We proceeded in the same way for the 29 
stars whose spectra did not reach beyond 2.4 pm. Now we 
could take advantage of all those stars which we extrapo¬ 
lated from 4 to 5 pm before in order to enlarge our input 
sample of stars for performing the interpolation. Like this, 
we ended up with spectra extending up to 5 pm for all IRTF 
stars. 

In the first part of Table |3] we present the spectral properties 
of the 180 spectra of our hnal stellar library which we will 
subsequently use in our SSP models (see next Section). 


8.1 Main ingredients 

W e calculated o u r SSP models based on the isochrones 
of iGirardi et al.l 1 I 2 OOOI) (hereafter: Pad o vaOO-i sochrones) 
and also the ones of IPietrinferni et all (l2004l l (BaSTI- 
isochrones). Both sets of model isochrones cover a wide 
range of ages and metallicities - in fact much wider than 
needed for our purposes given the limits in stellar atmo¬ 
spheric parameter coverage of the IRTF library (see Section 
n. Consequently, we calculated SSP models only for the 
hve metallicities [Fe/H] =—0.35,—0.25, 0.06, 0.15, 0.26 in 
the case of BaSTI and [Fe/H] = —0.40, 0, 0.22 for PadovaOO. 
Please note that the BaSTI isochrone corresponding to a 
metallicity of [Fe/H] = 0.15 is not yet publicly available since 
it was calculated e xclusively for this w ork (and others) by 
S. Cassisi (see also I Vazdekis et aLll2015l l. 

We also set some limits on the age range for which we ob¬ 
tained SSP models. So, we did not determine SSP models 
for ages younger than 1 Gyr, since at these ages, the con¬ 
tribution of young hot stars of Tea > 7000 K becomes im¬ 
portant. Unfortunately, the IRTF library is in general poor 
in these stars. The same holds for AGB stars which domi¬ 
nate the IR light at ages between 1 and 2 Gyr. In addition, 
this stellar evolutionary phase is difficult to model. Subse¬ 
quently, predictions differ a lot from each other depending 
on the synthetic prescriptions used by the aut hors. This can 
be see n by comparin g, e.g., the isochrone set of iMarigo et al.l 
ll2008ll to the one of lGirardi et al.l ll2000h . Apart from that, 
the IRTF library is also in general poor in luminous giants 
(AGB stars, see Section [S]), leading to a significant decrease 
in quality for SSP models of younger ages. 

We took various IMFs into account, obtaining models for 
IMFs of single slopes between 0.3 and 3.3. These slopes were 
either kept constant over the whole mass range or a shallower 
slope was adopted for masses smaller than 0.6 Mp) (so called 
unim odal and bimodal IMFs, respectively, IVazdekis et al.l 
Il996lf . In bimodal IMFs, the contribution of low mass stars is 
significantly reduced. We also calculated SSP models ba sed 
on the Kroupa and revised-Kroupa IMFs dKrou iiliool). In 
order to keep the following discussion more clear, we concen¬ 
trate on the SSP models and their predictions based on the 
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Kroupa IMF. This is justified since the infiuence of the IMF 
on our models is almost uegligible. Yet, possible impacts of 
different IMFs will be discussed in Section fl 1.2 1 


8.2 Computation of the models 

We present here a n extension of our models llVazdekis et al.l 
l2003l . I2OI0I . l2012fl . which are so far available for the opti¬ 
cal wavelength range between 3465 A and 9469 A, to the IR, 
out to 5/rm. The general modelling procedure follows very 
closely the approach carried out before in the optical. The 
general idea of this approach is to populate isochrones of 
various ages and metallicities with the spectra of our input 
stellar library (see Section [2]) according to the prescription 
given by a chosen IMF. Consequently, if the isochrone de¬ 
mands a star of a particular set of stellar parameters Teff, 
log((;) and [Fe/H], we try to reproduce this star as accu¬ 
rately as possible. For this, we use a so-called interpolator 
which creates a star of the desired parameters based on the 
input stellar library. 

The interpolator is a code which calculates a stellar spec¬ 
trum for any given set of stellar atmospheric parameters, i.e., 
effective temperature, surface gravity and metallicity, based 
on the stellar spectra of a whole library of input stars. The 
general idea of the interpolator code is to search for stars 
which are contained within up to eight cubes in the stellar 
atmospheric parameter space. One corner of these cubes is 
always centred on the parameters of the requested star. De¬ 
pending on how densely the parameter space is populated 
with stars in a certain region, the size of these cubes varies 
between a minimum and a maximum value. We adjust these 
input values of the interpolator to the scale required by the 
IRTF library. The parameters temperature and surface grav¬ 
ity are always taken into account, however, the metallicity 
only for stars of temperatures between 3250 and 7000 K. 
For hotter and cooler stars, respectively, we neglect this pa¬ 
rameter because the IRTF library does not contain a large 
enough number of stars in these temperature ranges in or¬ 
der to cover the parameter space in metallicity in a sufficient 
way (see Figure [5] and Section [5^ . 

After calculating the stellar spectrum of the requested pa¬ 
rameters, each of these spectra is normalised by convolution 
with the K band filter and the phot ometric librar ie s and 
colou r-temperature relations of, e.g., lAlonso et al. I (Il996l. 
Il999l j are used to convert the theoretical parameters (Teff, 
log((;), [Fe/H]) to the observational plane. We chose to nor¬ 
malise our stellar input spectra to the A'-band, since this 
band is located well in the middle of the spectral range cov¬ 
ered by the IRTF-library and spectrally close to the IR wave¬ 
length range between 2.5 and 5 /rm on which we focus in this 
paper. Finally, the stellar spectra are integrated along the 
isochrones. 

The procedure of calculating a SSP spectral energy distri¬ 
bution (SED) is summarised as: 

rmt 

Sa( t, [Fe/H]) = / Sxim,t, [Fe/H])• 

imi (3) 

A’(IMF, m, t) ■ FK{m, t, [Fe/H])dm, 

where SA(m, t, [Fe/H]) is the spectrum of a single star of 
mass m and a certain lifetime and metallicity which are 
determined by the adopted isochrone. Af(IMF,m,t) is the 


number of this kind of star which depends on the chosen 
IMF. The smallest and largest masses which exist for a sin¬ 
gle stellar population of a given age are designated by mi 
and mt- TW(m,t, [Fe/H]) describes the flux of each single 
stellar spectrum in the K band. We stress again that the 
fluxes obtained for the single stars requested by the theoret¬ 
ical isochrones arise from transformations based on empiri¬ 
cal photometric libraries and relations and not from the use 
of theoretical stellar atmospheres. This is one key difference 
compared to other models which make use of theoretical 
stellar atmospheres (see Section [TOl) . 


8.3 Quality and range of applicability 

In this section, we present a way to quantify the quality and 
reliability of our SSP models as a function of the parame¬ 
ters age and metallicity. To do this, w e adopt the normalised 
quality parameter Qn introduced in IVazdekis et al.l ll2010l l 
for the model predictions in the optical spectral range. The 
quality of an SSP depends on the parameter coverage of 
the input stellar library to be used by the interpolator (see 
Section EB in order to reproduce the single stellar spec¬ 
tra which are later integrated along the isochrone. A large 
number of stars with parameters differing as little as possi¬ 
ble from the ones requested for the computation of one re¬ 
spective model assures a high quality of the resulting model 
spectrum. The resulting quality value Q for one SSP model 
is defined as follows: 


Q = ' 



^7 = 1 "t 

NiFK^ 


c^oi y y 

/ «• \ 2 
[Fe/H]? ) 


j \^^IogC9)^ j 




Efci* NiFKi 


(4) 


In Eq. m Ni designates the total number of stars in the dif¬ 
ferent mass bins ranging from the smallest mass mi to the 
most massive stars of mass mt which can be encountered at 
the specific age t of the stellar population. In contrast. A/ 
stands for the number of stars which are found in each of the 
eight cubes in parameter space (see Section [6. HI . The sizes 
of these cubes are defined as and re¬ 
spectively, whereas ae ^, o-iog{g)m represent the 

minimum uncertainties of these three stellar atmospheric 
parameters. The minimum cube sizes are obtained by 
multiplying these a values by Xa. TWi gives the flux in the 
K band. Q is then normalised by dividing it by a minimum 
value Qm corresponding to the poorest still acceptable 
value of Q. For Qm, we assume that stars are located in at 
least three out of the eight cubes and that the maximum 
cube sizes do not exceed one-tenth of the total parameter 
ranges of our library. Consequently, models of values of 
larger than 1 are assumed to be of sufficient quality to be 
used safely. The interested reader can find a more detailed 
description of the definition of this parameter estima ting 
the quality of our SSP models in IVazdekis et all ll2010lj . 
Figure [5] displays the behaviour of the normalised quality 
parameter as a function of age for SSP models of differ¬ 
ent metallicities. In the case of underlying BaSTI isochrones 
(left panel of Figure O, the models of Solar and slightly sub¬ 
solar metallicities show the highest values of Qn exceeding 
the limit of 1 at least for ages larger than « 2 Gyr. For 
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Figure 8. Normalised quality parameter as a function of age 
for SSP models of different metallicities. The horizontal line at 
Qn = 1 marks the general limit of reliability of the models (see 
text). The left hand panel is dedicated to models based on BaSTI, 
whereas the right hand panel depicts the quality behaviour of 
PadovaOO based SSP models. 


PadovaOO based models (right-hand panel), only the SSP 
models of Solar metallicities provide valnes of Qn > 1 for 
the whole age range. The qnality parameters of the models 
of [Fe/H] = —0.40 already lie slightly below 1. Nevertheless, 
they can be considered acceptable as well. For lower metal¬ 
licities and also for supersolar valnes the models are of lower 
qnality and therefore must be handled with some care. This 
is equally true for the youngest ages between 1 and 2 Gyr. 
The low quality of the models in this age range can be also 
understood by the mismatch between the bolometric flux 
assigned to AGB and carbon stars in this age range and the 
actually available number of these stars in our input sample. 
Table [3] summarises the ranges of the various parameters for 
which we were able to compute reliable SSP models. Note 
that the values obtained here for the quality parameter Qn 
are much lower than the ones calculated for the model pre¬ 
dictio ns based on MILES (see Figure 6 in IVazdekis et all 
I 2 OIOI) . 


9 COLOURS MEASURED ON THE SSP SEDS 

In this section, we discuss the behaviour of colours measured 
from our SSP-SEDs as a function of age and metallicity. In 
order to obtain colours from our SSP model spectra, we 
convolve the spectral flux with the respective filter response 
functions. In this procedure, we calibrate the resulting mag¬ 
nitude s based on the Vega-spectrum of iGastelli fc Kurucd 
lll994h as a zeropoint. We calculated the magnitudes in the 
2MASS J, H , and Ks filters, in the Johnson J and K and in 
the Spitzer [3.6] and [4.5] /rm filters. As described in Section 
[T] we prepared all of the stellar spectra used in the modelling 
process to end at 5 p.m. However, the Spitzer [4.5] pm band 
extends slightly beyond 5 pm, up to 5.22 pm. Consequently, 
we need to quantify the missing flux in the [4.5] pm band. 
Therefore, we made use of a number of IRTF stars whose 
spectra extend well beyond 5 pm and calculated the flux in 
the complete Spitzer [4.5] pm filter as well as in a shorter 



Figure 9. Comparison between the behaviour with age of the 
{J—K) colour extracted from our SSP models and the same colour 
predicted from the m odels based on the MILES stellar library 
Uvazdekis et al . 1201 ( 1 ) . The dashed lines delineate the MILES pho¬ 
tometric predictions and the solid ones the colours calculated from 
our SSP models. The different colours represent different metal- 
licites (see legend). In the lower panel, we plot the absolute differ¬ 
ences between the colours. The dashed line marks the maximum 
difference of A = 0.02 mag which we consider still acceptable. 


]4.5] pm filter ending at 5 pm. The resulting difference in 
flux was only about A([3.6] — [4.5]) = 0.0001 and thus neg¬ 
ligible. So, the Spitzer ([3.6] — [4.5]) colour determined from 
our SSP models can be safely compared to the same colour 
from observations. 


9.1 Internal consistency of the model colors 

In Figure (9] we compare the (J — K) colour as calcu¬ 
lated from our model SSPs with the same colour origi¬ 
nating from the so-called MILES photometric predictions. 
These photometric colour predictions in the NIR have been 
determined by integrating the stellar fluxes obtained from 
the metallicity-dependent colour-temperature relations for 
dwarfs and giants described in Sections 15.21 and K2[ which 
are based on extensive empirical photometric stellar li¬ 
braries. Note that consequently, these colours were not ob¬ 
tained by convolving any synthesized SSP model SEDs with 
the respective filter res ponses, since the MI LES models ex¬ 
tend only up to 9500 A (IVazdekis et al.l[2012l ). Instead, these 
colours were calculated by using the metallicity-dependent 
colour-temperature relations mentioned above. 

As can be seen from Figure (H the (J — K) colour from the 
MILES photometric predictions agrees very well with the 
same colour measured from our SSP models over the whole 
range of ages, from 1 to 17.78 Gyr. Since these photometric 
predictions have only been published based on PadovaOO (see 
http://miles.iac.es), we present only comparisons to models 
using this set of isochrones. As is explicitly shown in the 
lower panel of Figure [9] the colours from our photometric 
predictions are less than 0.02 mag bluer than the colours in¬ 
tegrated from our SSP spectra for all ages larger than 3 Gyr. 
This holds for all metallicities within the covered range. For 
younger ages between 1 and 3 Gyr, the colour differences are 
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Table 3. Properties of the model SSP SEDs 


Spectral coverage 

Spectral range 

Spectral resolution 

Pixel size 

0.815 - 5 Atm 

(T = 60kms“^, (FWHM varies between 5.7 A at A = 0.945 Atm and 16.1 A at A = 3.475 Atm) 

2.73 A ■ pixel~^ 

Parameter coverage 

IMF type 

IMF-slope (for unimodal and bimodal) 

Isochrones 

metallicity (BaSTI) 

metallicity (PadovaOO) 

ages (BaSTI) 

ages (PadovaOO) 

unimodal, bimodal, Kroupa, Kroupa revised 

0.3-3.3 

BaSTI fPietrinferni et al. 2004'). PadovaOO fGirardi et al. 2000') 

—0.66dex < [Fe/H] < 0.26 dex 
-0.71 dex < [Fe/H] < 0.22 dex 

1 Gyr < t < 14 Gyr 

1 Gyr < t < 17.78 Gyr 


slightly larger. This close agreement shows that our models 
are internally consistent. 

9.2 Sensitivity of the Spitzer-^.Q\ — [4.5] colour to 
age and metallicity 

In this Section, we study the behaviour of the Spitzer 
([3.6] — [4.5]) colour as a function of age and metallicity. 
From Figure [TOl we see that this colour does not change much 
for ages above 2 Gyr. For a hxed metallicity, the Spitzer 
[3.6] — [4.5] colour becomes redder with increasing age by 
around 0.01 mag which is still comparable to the span of 
colour variation of around 0.002 mag cited by iMeidt et al.l 
(l2014li over the whole age range. This is equally true for the 
SSP models based on BaSTI as well as for the ones based 
on PadovaOO. For ages smaller than 2 Gyr, we observe a 
steep rise to redder colours. This is equally true for both the 
BaSTI and the PadovaOO based models, however it is more 
pronounced reaching redder colours in the former case. This 
increase in the Spitzer ([3.6] — [4.5]) colour can be explained 
by the enhanced contribution of AGB stars in this age range 
(see also fractions in Section 16.21) . The observed difference 
between the behaviour of the colours originating from the 
BaSTI and from the PadovaOO based models shows how sen¬ 
sitive this age regime is to the modelling of the AGB stars. In 
addition, the quality parameter (see Section IQ) decreases 
to values smaller than 1 for SSP models of ages younger than 
2 Gyr. 

Figure [in] also shows that the metallicity dependence of 
the Spitzer ([3.6] — [4.5]) colour is rather weak, however 
clearly visible. Over the whole age range, the overall dif¬ 
ference between the colours based on the isochrones of dif¬ 
ferent metallicities equals up to 0.04. Contrary to the be¬ 
haviour of the (J — K) colour with metallicity (see Figure 
O, the more metal-poor SSP models give slightly redder 
Spitzer ([3.6] — [4.5]) colours compared to the more metal- 
rich ones. A likely explanation for this reversed trend is the 
presence of the CO absorption band in the [4.5] /rm Spitzer 
filter (see also Sectio n |3. This obse r ved tr end i s also seen by 
[Peletier et al.l (l2012ll : iMeidt et al.l (l2014l l and I Norris et al.] 
(l2014ll . who relate it to the absorption by the CO in the 
[4.5] pm band. Larger metallicities lead to increased CO ab¬ 
sorption, and subsequently to a smaller flux in the [4.5] pm 
band. This, however, makes the Spitzer ([3.6] — [4.5]) colour 
slightly bluer compared to the same colour measured from 


SSP models of lower metallicities and thus with less CO 
absorption. This trend just does not hold for the mod¬ 
els of [Fe/H] = —0.35 from which we measure a Spitzer 
([3.6] — [4.5]) colour which is slightly bluer than the one 
determined from the models of [Fe/H] = —0.25. The differ¬ 
ence however is small (<= 0.01 mag). Since also the quality 
of the models of a metallicity of [Fe/H] < —0.3 mag is only 
modest (see Figure [S] in Section [8.311 . we are not in a posi¬ 
tion to decide whether we observe here a real saturation of 
the Spitzer ([3.6] — [4.5]) colour with decreasing metallicity 
or not. Our data do not permit us to investigate in detail 
the depth of the CO band as a function of metallicity, as 
the IRTF spectra show a gap due to telluric absorption in 
this range, which was corrected artificially using the theo¬ 
retical Phoenix library (see Section [7||. Yet a comparison to 
the colours obtained from other (theoretical) models in this 
wavelength range which reach lower metallicitie s seems to in¬ 
dicat e a continuation of this trend of reddening llNorris et al.l 
|2014 , see below). 

The small metallicity dependence of the Spitzer ([3.6] — [4.5]) 
colour is in agreement with the literature. iMeidt et al.l 
ll2014h give a variation of A([3.6] — [4.5]) = 0.06 over the 
larger metallicity range from [Fe/H] = —1.7 to [Fe/H] = 0.4. 
This corresponds to A([3.6] — [4.5]) = 0.01 over the metallic¬ 
ity range A [Fe/H] = 0.35, which is the ran ge of highest reli¬ 
ability of our models. IPeletier et al.l ll2012|j state that barely 
any fluctuation of the Spitzer ([3.6] — [4.5]) colour is vis¬ 
ible for metallicities higher than Z = 0.008, which equals 
[Fe/H] = -0.4. 

Finally, we integrate the spectral flux from our models falling 
within the two filters W1 and W2 o f the Wide-held Inf rared 
Survey Explorer (WISE) satellite ilWright et al.ll20f^ and 
thus obtained the WISE {W1 - W2) colour. The WISE W1 
and W2 filters are only slightly shifted with respect to the 
Spitzer [3.6] pm and [4.5] pm bands. Subsequently, studying 
the behaviour of the WISE {W 1 — IT2) colour with age (see 
Figure [TT} allows us to confirm the trends observed for the 
Spitzer colour (see Figure [Toll . Indeed, when comparing Fig¬ 
ure [To] to Figure [TT] we can conclude that these two colours 
do behave basically the same as a function of age and metal¬ 
licity. Due to the slightly different filter responses which are 
shifted with respect to each other, the WISE (W1 — W2) 
colour is generally redder by about 0.03 mag for ages younger 
than 2 Gyr and by still about 0.01 — 0.02 mag for older ages. 
Furthermore, the WISE W2 filter seems to be slightly less 
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Age [Gyr] Age [Gyr] 


Figure 10. Behaviour of the Spitzer ([3.6] — [4.5]) colour deter¬ 
mined from our SSP models with age. The different colours and 
symbols represent different metallicites (see legend). The colour 
predictions corresponding to the highest and lowest metallicity 
values are of slightly lower quality. The left panel displays colours 
from models based on BaSTI, whereas the colours in the right 
one originate from models based on PadovaOO. 



Figure 11. Same as Figure ITol but for the WISE (W1 — W2) 
colour. 


affected by the CO absorption band centred at 4.6 fim, since 
the blueing of colours with increasing metallicity is a bit 
less pronounced. In a recent publication. iNorris et al.l (|2014l i 
have determined WISE photometry of a diverse sample of 
globular clusters and early-type galaxies and then compared 
the resulting (W1 — W2) colour to the predictions of some 
stellar population models (see next Section [iQl). As can be 
seen from their Figure 4, their observed sample shows the 
same blueing of the WISE {W1 — W2) colour with increasing 
metallicity which we observe from our models. Apparently, 
this trend holds even for lower metallicities which we can¬ 
not reach with our models due to the limitation of the IRTF 
library. 


10 COMPARING TO OTHER SSP MODELS 
FROM THE LITERATURE 


In this Section we compare our predictions with 
other models ava il able in the literature. According to 
iBarmbv fc JalilianI ll 2012 li . there exist only three other 
models which provide explicit colour predictions in the 
wavelength range beyond the K band. These three are 
the SSP models of the Padova group emer ging as a by¬ 
prod u ct of their library of stellar isochrones ( Marigo et al.l 
I2OO8I : iGirardi et al.l 2 010 l. here after Marigo-models), the 
GALEV mod els (|Kotulla et al.l l2009l l and the FSPS ones 
jConrov et al.l[2009l . 2010h . However, it should be mentioned 
that IR colours can be also obtained from all other models 
which are based on the theoreti c al BaSeL library (iKurucd 
Il99^: iLe ieune et all Il99 7l. Il998l : IWestera et al.l l^02l ~) like 


Bruzual fc Chariot! ll2003[ l and Maraston ( 20051 ) by con- 


e.g. 

volving the predicted spectra with the respective filter re¬ 
sponse functions. In what follows we focus however on the 
three models g iving directly colour predic tions in the IR dis¬ 
cussed also bv iBarmbv fc JalilianI (l2012t) complemented by 
th e predictions o f the n ewest release of the Padova models 
bv iBressan et al.l (120121 '). 

All of these models are based on Padova isochrones. How¬ 
ever, the models differ from each other in the stellar li¬ 
braries on which they are based. Contrary to our mod¬ 
els, all of the underlying stellar libraries are theoretical 
ones. So, the GALEV model s use version 2.2 of the theo¬ 
retical BaSeL stellar library ilLeieune et al.l[l998h . whereas 
the FSPS models rely on the updated vers ion 3.1 of this 
same spectral library JWestera et al.l I 2 OO 2 I I in the wave¬ 
length range between 2.5 /rm and 5 /rm. The spectr a of the 
BaSeL library are based on stellar atmospheres fromiKurucz 

1992) and - in the case of M giant s - from iFluks et al. 

1994 ) and iBessell et al.l d 19891 . Il99ll) . The Marigo mod¬ 


els main ly make use of th e revis ed ATLAS9 model spec¬ 
tra from ICastelli fc Kurucd ll2003l) . For M, L and T dwarfs, 
they adopt the model at mospheres o f Allard et al.l (l2000l) 
and for M giants those of lFluks et al.l (Il994l) . Moreover, the 


prescriptions for the thermally pulsating asymptotic giant 
branch (TP-AGB) stars, vary between the models. In or- 
der to account for th e m, FS PS uses empirical spectra from 
iLancon fc Mouhcind ll2002h and the Marigo models syn¬ 
thetic spectra from Loidl et (l200ll) . As these two spectral 
libraries extend only up to [2.5] /rm, they have to be ex¬ 
trapolated with a Rayleigh-Jeans model. For carbon stars, 
the FS PS models are based on the models bv lAringer et al.l 
(l2009l) . The GALEV models do not treat post-main se¬ 
quence stars in a particular way. For all of these four fami¬ 
lies of models, we obtained the Spitzer [3.6] fim and [4.5] pm 
magnitudes for a Solar metallicity and for all available ages 
older than 1 Gyr. The models adopted by us are based on a 
Kroupa-like IMF. 

In the case of the GALEV models, we used the default 
parameters and selected to neglect the dust contribution. 
For the Marigo mod els, we used the e volutionary tracks by 
iMarigo et al.l (l2008l ) together with the IGirardi et all (|201(]|) 
case A correction for low-mass, low-metallicity AGB tracks, 
whereas we used ve rsion 1.1 of the PARSEC models of 
iBressan et al.l ll2012l ). For both sets of models, we did not 
select any circumstellar dust. In the case of FSPS, we re- 
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trieved the models made readily available online based on 
Padova isochrones and the BaSeL stellar library. 

Figure [12] shows the behaviour of all these models together 
with our new SSP models based on BaSTI and PadovaOO 
isochrones, respectively, for ages older than 1 Gyr. The 
Spitzer ([3.6] — [4.5]) colours predicted by the Marigo mod¬ 
els and the Bressan models agree almost perfectly with our 
two new sets of models over the whole age range. For ages 
older than about 2 Gyr, the Marigo and the Bressan models 
give colors which are generally only around 0.005 mag redder 
than those of our two sets of models. 

The Spitzer ([3.6] — [4.5]) colours derived from the FSPS 
models also coincide reasonably well with the colours de¬ 
termined from our sets of SSP models and the Marigo and 
Bressan ones. The almost constant colours obtained from 
FSPS are about 0.01 — 0.015 mag redder than the ones mea¬ 
sured from the former ones for ages older than 2 Gyr. For 
the youngest ages, the FSPS models do not show the turn 
to redder colours due to the enhanced contribution of AGB 
stars which become s visible in our models and in the Padova 
ones (compare also lMarastonlflO^ . 120051 '). 

The GALEV models, however, give absolute values of the 
Spitzer ([3.6] — [4.5]) colour which are completely offset to 
redder colors by about 0.10 — 0.13 from our SSP models. 
Most likely, there are two main reasons why this set of mod¬ 
els fails to provide colour predictions which lie in the same 
range as all the others, although their input parameters do 
not differ much from e.g. the Marig o or Bressan models (see 
above and iBarmbv fc Jalilianl[2012l l. On the one hand, their 
far-too-red colours might be explained by their choice of the¬ 
oretical stellar atmospheres for transforming the theoretical 
stellar parameters to the observational plane. On the other 
hand, they do not take the GO absorption in the Spitzer 
[4.5] fim band into account. The crucial point of correcting 
for this absorption in order to obtain re liable colour predic - 
tions in the IR is also pointed out by I Norris et al.l ll2014l l 
who compare their observations also to the GALEV mod¬ 
els (see_tliei£^]igure_^.Th^ claim that the FSPS models 
of IConrov et aLl ll2009l . l20ld l do not consider the GO ab¬ 
sorption neither, a finding that we cannot confirm from the 
predictions visuali sed in our Figure 1121 Ap art from that, 
I Norris et al.l (l2014h state that the models bv iBressan et al] 
(I 2 OI 2 II are the first modern SSP models which are able to 
reproduce the observed {W1 — W2) colours of dust-free stel¬ 
lar populations. Since the colour prediction s arising from 
our SS P models closely agree with the ones of iBressan et al.l 
(l2012tl . we can argue that this also holds for our newly cal¬ 
culated models based on empirical stellar spectra. 


11 MASS-TO-LIGHT RATIOS IN THE 

SPITZER [3.6] AND [4.5] BANDS AND 
DEPENDENCE ON THE IMF 

In this Section, we study the behaviour of the mass-to-light 
ratios in the Spitzer [3.6] pra and [4.5] /rm bands as a func¬ 
tion of age and metallicity. We also investigate the impact 
of varying the IMF on the mass-to-light ratios in those two 
bands. As expected, here the effect of the IMF is much more 
pronounced than for the colours. Nevertheless, we also dis¬ 
cuss briefly the IMF dependence of the Spitzer ([3.6] — [4.5]) 
colour (see Subsection lll.2l) . 



Figure 12. Comparison of the behaviour of the Spitzer ([3.6] — 
[4.5]) colour with age for four different models (see legend) and 
our own (BaSTI- and PadovaOO based) ones. All of the displayed 
models are of Solar metallicity. 


11.1 Mass-to-light ratios in the IR 


Apart from the spectra of the SSPs, our models also give 
their stellar masses as an output. Together with the flux in 
a specific band which can be obtained by integrating over the 
model spectra, we are subsequently able to calculate mass- 
to-light ratios in various bands. In what follows, we deter¬ 
mine the mass-to-light ratios in the Johnson K band as well 
as in the Spitzer [3.6] pm and [4.5] pm bands. We provide all 
of our mass-to-light ratios in Solar units, Mq /L©. Therefore, 
we adopted the absolut e Solar magnitud es of Mq® = 3.24 
and M|i'® = 3.27 from Oh et al.l feoosll and Mq = 3.28 
from iBinnev &: Merrifieldl lll99^ in order to calculate the 
luminosities. 

The mass-to-light ratios measured in these NIR and IR 
bands present the great advantage that they depend less on 
parameters like age and metallicity and hence are more sta¬ 
ble and constant comp ared to their counterpa rts determined 
in the optical regime. iBell fc de Jond (l200l[l . e.g., mention 
that over the same parameter range, the mass-to-light ra¬ 
tios change by a factor of about 7 in the B band, 3 in the 
I band and 2 in the K band. This behaviour of the mass- 
to-light ratio can be understood considering that the light 
of SSPs in the IR is dominated by red giant branch (RGB) 
stars. With advancing time, the RGB gets more and more 
populated by stars of later spectral type with decreasing 
mass-to-light ratios. At the same time, however, the mass- 
to-light ratio of the remaining main sequence stars increases 
since it is mainly determined by later type M stars of higher 
mass-to-light ratio. Consequently, the overall changes in the 
mass-to-light ratio with time are rather smooth as these tw o 
effects compensate each other (see, e.g.. lMeidt et al.]l2014fl . 
Such a generally small parameter dependence is also con¬ 
firmed by Figure [T51 depicting the M/L ratios in the three 
NIR and IR bands and for three (in the case of SSP mod¬ 
els based on BaSTI isochrones) and two (SSP models based 
on Padova isochrones) different metallicities. In particular, 
there is basically no dependence on met allicity. This insensi - 
tivity to metallicity is also confirmed bv iNorris et al.l (l2014j i 
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Figure 13. Behaviour of various M/L ratios with age. Blue lines 
represent the M/L ratio in the K band, black lines in the [3.6] pm 
band and red lines at [4.5] pm. The different line styles stand 
for the metallicities as indicated in the legend. The M/L ratios 
plotted in the left hand panel originate from SSP models based 
on BaSTI isochrones, whereas the ones shown in the right hand 
panel were obtained from SSP models based on PadovaOO ones. 


who calculate M/L ratios in the WISE W1 and W2 bands 
from the newest version of the Padova models llBressan et al.l 
[2 oT 2). The differences between the M/L ratios in models of 
varying metallicity are smaller than the differences between 
the M/L ratios measured in the various filters. However, 
the dependence with age is smallest for the M/L ratio in 
the [3.6] pm band, whereas the time evolution of the M/L 
ratio in the [4.5] pm band resembles a lot the one in the 
K band (see again Figure [T3|. The M/L ratios of the two 
latter ones are slightly larger than the one measured in the 
[3.6] pm band. In general, the M/L ratio determined in the 
[3.6] pm Spitzer band just changes from {M/L) 3 .q « 0.20 
for an age of 1 Gyr to {M/Ljs.e « 1.1 at the maximum age 
of 17.78 Gyr, i.e., by a factor of 5.5. iNorris et al.l ll2014l ') find 
an increase of the M /L-ratio in the WISE bands by a factor 
of 2 between 3 and 10 Gyr which agrees perfectly with our 
results (compare Figure 031 to their Figure 5). In contrast to 
this, the M/ L ratio in the V ban d calculated from MILES 
SSP models (IVazdekis et al.ir2010l l of Solar metallicity grows 
from {M/L)v ~ 0.6 to {M/L)v ~ 6 within the same age 
range, i.e., it increases roughly by a factor of 10. 
Furthermore, comparing the left and right hand panels of 
Figure [THl it can be seen that theM/L ratio mildly depends 
on the employed isochrones. 

On the contrary, the M/L ratios are highly sensitive to the 
IMF used in the modelling procedure, as Figure [U clearly 
depicts. In general, M/L ratios originating from SSP mod¬ 
els based on unimodal IMFs are larger than the ones based 
on bimodal IMFs. This is due to the enhanced contribu¬ 
tion of faint low-mass stars in the case of unimodal IMFs 
compared to bimodal IMFs for which the number of stars 
wi th masses smaller than 0.6 Mq is reduced (e.g., Figure 2 
in iLa Barbera et ah I [113). Also note that a bimodal IMF 
of slope 1.3 is virtually identical to the Kroupa IMF. More¬ 
over, it can be seen from Figure [2] that a steepening in 
IMF slope leads to a greater increase in the M/L ratio 
for unimodal IMFs than for bimodal ones. According to 



Figure 14. Behaviour of various M/L ratios at [3.6] pm with 
age. The depicted M/L ratios originate from SSP models of So¬ 
lar metallicities calculated using different IMFs (see legend). The 
M/L ratios plotted in the left-hand panel are calculated from 
SSP models based on BaSTI isochrones, whereas the ones shown 
in the right-hand panel were obtained from SSP models based on 
PadovaOO isochrones. 


iLa Barbera et al] ll2013l) . the M/L ratios arising from mod¬ 
els based on unimodal IMFs give too large values, exceed¬ 
ing by far dynamically determ ined M/L ratios like, e.g., 
those of iGappellari et al] (l2013l ). Stellar M/L ratios calcu¬ 
lated based on bimodal IMFs, however, are consistent with 
the latter ones also in the most bottom-heavy c ases (i.e., for 
the largest s l opes, see also ISpiniello et ^l2014l ). 
iMeidt et al.l (l2014l ) c alculated (M/Lja.e r atios based on the 
model predictions of iBruzual fc Gharlotl (l2003l ) which they 
then transformed to ([3.6] — [4.5]) using a newly calibrated re¬ 
lation between the NIR and the IR colors of giants observed 
with GLIMPSE l|Beniamin et al.l l200,3l : IChurchwell et al] 
l2009t) . They find only very small variations of the (M/L/s.e 
ratios as a function of age and metallicity. That is why they 
suggest to use even a constant M/L ratio of M/L^.q ~ 0.6 
independent of the Spitzer ([ 3.6[—[4.5]) colour and conse- 
quently of age and metallicity (|Meidt et al.ll20l3 ). Recently, 
iNorris et al. (|2014l ) confirmed that based on this constant 
M/L ratio, it is possible to determine stellar masses with 
errors no larger than 0.1 dex in the middle of the age range 
3 to 10 Gyr. The value of M/L 3.6 « 0.6 is also quite well 
compatible with an average value for our M/L ratios. How¬ 
ever, we obtain a larger dependency on the two parameters 
age and metallicity, which is also due to the fact that our 
models are single-burst like and do not contain expone ntially 
declining star formation histories l|Norris et al.ll20141 . 

11.2 IMF dependence of our models 

Up to here, we have always discussed the SSP models based 
on a Kroupa-like IMF. However, as we pointed out in Section 
18.11 we also calculated our SSP models for a suite of IMF 
shapes and slopes. The two sets of IMFs which we used 
were the unimodal and the bimodal ones fsee lVazdekis et ^ 
Il996h . The varying behaviour of a unimodal and bimodal 
IMF of the same slope p for stellar masses smaller than 






























































16 B. Rock et al. 


aroun d 0.5 Mq is visualized in Figure 3 of iLa Barbera et al.l 
ll2013li . 

The two upper panels of Figure [15] show the effect of vary¬ 
ing the IMF slope in the case of bimodal IMFs. The corre¬ 
sponding figure for unimodal IMFs is not shown, but looks 
essentially the same. The general trend that becomes visi¬ 
ble is a slight reddening of the Spitzer ([3.6] — [4.5]) colour 
with increasing IMF slope of the order of 0.02 — 0.03. How¬ 
ever, the difference in colour between two models differing 
by A(]Fe/H]) « 0.40 is comparable to the difference between 
two models of the same metallicity but a A{fi) of 1 — 1.5. 
Consequently, we can conclude that the effect of varying the 
IMF slope p. on the colours is small for our models. 

In the two lower panels of Figure [15] we compare the colour 
predictions from SSP models based on a bimodal and a uni¬ 
modal IMF, respectively, for the same slopes. We find that 
a bimodal and a unimodal IMF of slope 1.3 give virtually 
the same Spitzer ([3.6] — [4.5]) colour over the whole age 
range. In the case of the steeper slope of 2.3, a small differ¬ 
ence of again 0.02 — 0.03 arises. The colours deduced from 
SSP models based on unimodal IMFs become increasingly 
redder compared to the ones calculated from models with 
underlying bimodal IMFs. However, the observed effect is 
again rather small. 

Finally, we can summarise that the impact of varying both 
the IMF type and slope on the behaviour of the Spitzer 
([3.6] — [4.5]) colour measured from our SSP models is rather 
small. The observed dependences described above are of the 
same order as those seen when varying age and metallic¬ 
ity for ages older than ~ 2 Gyr and over our limited range 
of metallicity (see Section 0. In general, we can conclude 
that the Spitzer ([3.6] — [4.5]) colour can be used to a cer¬ 
tain extent as a proxy for distinguishing populations of ages 
younger than 2 Gyr from older ones and also more metal- 
rich populations from more metal-poor ones. However, it 
remains very difficult to make a clear statement about the 
underlying IMF of a stellar population based on this colour. 


12 COMPARISON TO OBSERVATIONS 

In this Section, we carry out a comparison between the 
colours determined from our SSP models and observed pho¬ 
tometry to validate our models. We try to reproduce the 
colours (see Section [12.Ill and the colour-velocity dispersion 
relations (see Section [12.211 of early-type galaxies. For this, 
we make use of the sample of 48 galaxies observed with the 
Spectrograp hic Areal Unit for Research on Optical N ebulae 
(SAURON) ([Bacon et al.ll2nnil : Ide Zeeuw et al.ll2nn2h . 

12.1 Early-type galaxies 

Figure [16] compares the observed Spitzer ([3.6] — [4.5]) 
colour of the 48 nearby, early-type elli ptical and lenticu¬ 
lar galaxies from the SAURON survey (iBacon et al.l l2(lf)ll : 
Ide Zeeuw et al.]]2002l l to the predicted colours of our SSP 
models as a function of age. The age here is the optical SSP 
age determined from th e SAURON data within 1/8 • reft by 
]Kuntschner et ahl (]201Cll l. These data show that galaxies be¬ 
come gradually bluer as they become more massive. For this 
plot, the Spitzer colours were measured inside 1/8-reft, ren¬ 
dering them bluer than when determined inside the whole 


effective radius reft. Both effects are consistent with metal¬ 
licity decreasing outwards in the galaxies, causing the CO 
band to weaken and the ([3.6] — [4.5]) colour to redden. Our 
models match quite well the observed colours of the galax¬ 
ies with the oldest mean luminosity-wei ghted ages, as de¬ 
term ined from the optical spectral range i Kuntschner et aP 
I 2 OIOII . However, they are unable to reproduce the redder 
colours of younger galaxies. These galaxies are likely to have 
more extended SFHs and thus cannot be modelled as a 
single stell ar population . This is confirmed by the classi¬ 
fication of ]Shapiro et al.1 (]201Clfl who estimate the level of 
star formation present in the SAURON galaxies based on 
the non-stellar dust emission in the [8] pm Spitzer band. 
The galaxies for which they report clear and likely star for¬ 
mation, respectively, coincide with the reddest galaxies in 
Figure 1161 Most of the non-stellar emission in the Spitzer 
[3.6] pm and [4.5] pm bands is due to poly cyclic aromatic 
hydrocarbons ('PAHsf. ISchodel et al.1 (1201411 studied the ra¬ 
tios of the diffuse emission, i.e., the remaining emission after 
point source subtraction, in the Spitzer [3.6] pm, [4.5] pm 
and [8] pm bands for observations of the Galactic centre. 
They estimate a PAH contribution of 0.03 times the flux at 
[8] pm in Jansky for the [3.6] pm band and of 0.04 times the 
[8] pm flux for the [4.5] pm band. We used these estimates 
in order to correct the Spitzer ([3.6] — [4.5]) colours of the 
SAURON galaxies for PAH emission. As can be seen from 
Figure [T71 we end up with slightly bluer colours which are in 
better accordance with our models. This is obviously only a 
crude estimate. Moreover, the colours displayed in Figure ITTI 
were measured within one effective radius contrary to Figure 

M 

Figure [18] displays the Spitzer colours of the same galaxies, 
but this time plotted versus metallicity, which was also de¬ 
termin ed from the SAURON spectra by iKuntschner et aD 
1 I 2 OIOII . Like in the previous Figure 1161 the colours and 
metallicities were measured within 1/8 • Teff. We overplot 
two models of ages 2 and 10 Gyr, based on both BaSTI and 
PadovaOO, respectively. Again, our SSP models coincide best 
with the bluest of the SAURON galaxies, whereas the red¬ 
der, younger, mostly star-forming galaxies lie outside the 
range traceable by our models. Figure [TS] hints that the ob¬ 
served redder ([3.6] — [4.5]) colours are rather caused by ex¬ 
tended SFHs including also younger populations (compare 
to Section than due to lower metallicities. After all, 
according to Figure 1181 the metallicities of the SAURON 
early-type galaxies are at most slightly subsolar. However, 
in this paper, we refrain from drawing a final conclusion 
concerning this point. 

12.2 Confirming the observed colour-velocity 
dispersion relations of galaxies 

]Ferreras et al.1 (]2013fl and ]La Barbera et ahl (]2013fl analyse 
a variety of optical IMF-sensitive spectral features in a 
large sample of early-type galaxies and draw conclusions 
on the underlying IMF. The former ones additionally carry 
out a spectral fitting using the MIU SCAT population syn¬ 
thesis models (IVazdekis et al.1 12 OI 2 II . Both find clear evi¬ 
dence for a relation between the IMF and the central ve¬ 
locity dispersion (jq. According to their studies, galaxies of 
(Jo ~ 120 kms“^ are characterised by a standard Kroupa-like 
IMF, whereas galaxies showing a central velocity dispersion 
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Figure 15. Behaviour of the Spitzer ([3.6] — [4.5]) colour as a function of age for SSP models of different IMF types and slopes. The 
upper left and right panels display the Spitzer ([3.6] — [4.5]) colour for bimodal IMFs of the three different slopes p, = 1.3, 2.3, 2.8 (see also 
legend). The blue lines show the predicted colours for a subsolar metallicity and the green lines for the Solar one, respectively. Models 
based on BaSTI are presented in the left-hand panel, whereas PadovaOO based models are shown in the right-hand one. The lower two 
panels visualise the Spitzer ([3.6] — [4.5]) colour originating from SSP models based on Solar metallicities and uni- and bimodal IMFs of 
the two slopes p = 1.3, 2.3. 


of ao ~ 200kms“^ and even larger, are best described by a 
bottom-heavy bimodal IMF with slope oi p > 2. This be - 
haviour is visualised i n Figure 12 ofiLa Barbera et al.l (l2013l) 
and i n Figure 4 of iFerreras et al1 ~i 2013l )~ Spiniello et ahl 
l|2014h hnd a similar correlation between the central velocity 
dispersion and the IMF slope for early-type galaxies. 


iFalcon-Barroso et al.1 l|2nill ) present tight relations between 
the {V — [3.6]) colour and the central stellar velocity dis¬ 
persion CTe measured within one effective radius for a sam¬ 
ple of elliptical and SO galaxies from the SAURON survey 
(see their Figure 7, left bottom panel). In the middle and 
bottom panels of Figure [1^ we re produce the same co lour 
versus cr plot for the sample from IPeletier etThl (|2012| . see 
their Figure 11a). Additionally, we overplot the predicted 
colours from our SSP models for two stellar populations of 
different cr. We assume a small a of 120kms“^ and a larger 
on e of 230kms~^ . resp ectively. According to the relation 
of ISpiniello et al.l (l201^ ). we choose for the model galaxy 
of (T = 120kms“^ a Kroupa IMF of /r « 1.3, whereas a 
galaxy of ct = 230 km s”*^ is described by a bimodal IMF 
oi p ~ 2.3. Furthermore, we assumed Solar metallicity and 


adopted an age of 6 Gyr for the former and of 12 Gyr for 
the latter galaxy. Our V band prediction comes from the 
MILES photometric predictions (see Section [3)l , since the V 
Hlter lies outside the wavelength range covered by our new 
SSP models. 

From Figure [19] (upper panel), we can deduce that the lo¬ 
cation of our model galaxies in the {V — [3.6]) — cr space is 
in very good agreement with the data. The same statement 
holds when plotting the (V — [4.5]) colour versus the veloc¬ 
ity dispersion cr, as seen in the middle panel of Figure 1191 
However, as seen in the lower panel of Figure [T^ we do not 
obtain such a good agreement between the observations and 
the modelled data points for the Spitzer ([3.6] — [4.5]) colour 
versus a. 

The reddening of the colours of the galaxies as cr decreases 
is not reproduced well by our SSP models, for which the 
predicted ([3.6] — [4.5]) colour is clearly too blue. This is ob¬ 
viously the same effect as observed before in Section 112.11 
There, we explained it by the extended SFHs of the younger, 
low-mass galaxies which most likely include a signihcant 
fraction of thermally pulsating AGB and other young stars. 
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Figure 16. Spitzer ([3.6] — [4.5]) colour versus age for 48 nearby 
early-t ype elliptical and lentic ular galaxies from the SAURON 
survey fauntschner et ^|20101) . Errorbars denote the uncertain¬ 
ties of both colour and age. Galaxies with clear and possible signs 
of ongoing star formation are shown as red triangles and green 
squares, respectively, whereas those without star formation are 
delineated as black dots. The coloured lines show the predictions 
from our SSP models for various metallicities. Solid lines corre¬ 
spond to models based on BaSTI, whereas the dashed ones delin¬ 
eate PadovaOO based models. 



Figure 17. Same as Figure fTbl but here age and Spitzer colour 
are determined inside one effective radius. The black dots de¬ 
lineate the SAURON galaxies before they were corrected for 
PAH emission. The PAH-corrected galaxies are displayed with 
red squares. 


This finding is in accordance with iFalcon-Barroso et al.l 
ii2oiiri . who report a large contribution of these kind of stars 
to the overall galaxy light in the IR. Contrary to the situ¬ 
ation in the optical wavelength range, these types of galax¬ 
ies apparently cannot be reproduced by a single-burst stel¬ 
lar population model in the IR. Consequently, this wave¬ 
length range might enable us for the first time to disentangle 
the various stellar populations present in young, low-mass, 
metal-poor galaxies as well as to estimate the contribution 
of the non-stellar PAH emission. Further investigation on 



Figure 18. Same as Figure [T6l but now plotted versus metallicity. 
Errorbars denote the uncertainties of both colour and metallicity. 
The blue and green lines depict SSP models of 2 and 10 Gyr, 
respectively. 



Figure 19. Colour versus a visualisations of a selected sam¬ 
ple of nearby early-type galaxies from the SAURON survey 
(red squares). The overplotted large blue dots are the predicted 
colours from our BaSTI based SSP models for bimodal IMFs of 
slopes p = 1.3 and p = 2.3, corresponding to the two values 
of cr = 130kms“^ and a = 230kms“^, respectively. The black 
straight lines delineate the best fits to the observed galaxies. The 
upper panel shows the (V —[3.6]) colour, the middle one (U—[4.5]) 
and the lower one the Spitzer ([3.6] — [4.5]) colour. 


the different stellar populations and non-stellar components 
in these galaxies will be studied in a future paper. 


13 SUMMARY AND CONCLUSIONS 

We present new evolutionary single stellar population mod¬ 
els covering the so-far poorly studied IR wavelength range 
beyond the K band, i.e., between 2.5/rm and 5 fj,m. Our 
models use as stellar input 180 carefully sel ected and pre¬ 
pared spectra from the I RTF stellar library (ICushing et al.l 
I 2 OO 5 I : iRavner et al.l[200^ . We have determined that the res¬ 
olution of the IRTF spectra is 60kms“^ (see Table |3]). This 
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makes them very appropriate to study most elliptical galax¬ 
ies. We determined the atmospheric parameters of the stars 
as accurat ely as possible using a compilation from the liter- 

and the colour-temperature 


ature (see [Cenajro_^£_^ 


relations of Alonso et al 


using 

1200 '^ 


( 1996l.ll999li . 


We provide the model spectra in a way which enables the 
reader to obtain all kinds of R colours by convolving them 
with the respective filter response functions. Since the stars 
of the IRTF library are mostly bright, nearby stars of around 
Solar metallicity, our SSP models are most accurate in the 
metallicity range between [Fe/H] = —0.40 and [Fe/H] = 0.2. 
We also obtained model predictions for slightly higher and 
lower metallicities, which are, however, of lower quality. 

We calculated our models based on both BaST I 
llPietrinferni et ^|2004^ and PadovaOO llGirardi et al.l[2005 l 
isochrones. They mainly differ for ages younger than around 
3 Gyr, since they use different prescriptions for the AGB star 
contri bution whi c h bec o mes c rucial at these ages. Accord¬ 
ing to iMarastonI (Il998l . l2005ll . AGB stars are particularly 
important at ages between 1 and 2 Gyr. The increasing con¬ 
tribution of this phase results in redder ([3.6] — [4.5]) colours. 
We cannot extend our models to ages younger than 1 Gyr 
due to an insufficient number of AGB stars and young, hot 
stars in the employed IRTF library. A summary of most 
of the points mentioned so far can be found in Table |3] in 
Section 18.31 

We obtained our SSP models for Kroupa-like, revised 
Kro upa and unimodal an d bimodal IMFs of varying slopes 
fsee IVazdekis et'^l2003l . for the different IMF types). The 
slopes of the latter two cases range from ^ — 0.3 to /r = 3.3. 
A unimodal IMF of slope 1.3 resembles a Salpeter-IMF and 
a bimodal one of the same slope a Kroupa-IMF. As our 
analysis (see Section 111.211 shows, our models basically do 
not depend on the used IMF. Gonsequently, all the other 
studies presented in this paper are carried out for models 
with an underlying Kroupa-like IMF. 

In what follows we summarise the most important results 
which we obtain from our models. 


(i) We tested the internal consistency of our model predic¬ 
tions by comparing the (J—K) colour measured on the IRTF 
based SSP spectra to the photometric predictions based on 
extensive photometric libraries published on our webpage 
(http://miles.iac.es) and found a good agreement, within 
0.02 mag over most of the age range. 

(ii) In general, we can observe only a weak dependence 
on age and metallicity for our SSP models. The Spitzer 
([3.6] — [4.5]) colour is apart from a very small trend of 
reddening almost constant with age for ages older than 
« 3 Gyr. In this age range, the small variations in the Spitzer 
([3.6] —[4.5]) colour are almost exclusively driven by metallic¬ 
ity. For younger ages between 1 and 3 Gyr, the ([3.6] — [4.5]) 
colour shows some gradient. It becomes significantly redder 
depending on the used isochrones (BaSTI or PadovaOO) due 
to the enhanced AGB star contribution. Solar and supersolar 
metallicities result in bluer ([3.6] — [4.5]) colours than lower, 
subsolar ones. This at hrst sight unexpected behaviour can 
be explained by the prominent CO absorption band in the 
[4.5] pm filter which becomes more important for stars of 
lower temperature and larger metallicity. Hence, for these 
metallicities, a greater amount of flux is absorbed compared 
to lower ones causing the observed bluer ([3.6] —[4.5]) colour. 


Moreover, the ([3.6] — [4.5]) colour is found to be almost un¬ 
affected by the underlying IMF. 

(iii) The ([3.6] — [4.5]) colours predicted by our SSP mod¬ 
els lie in the same range as the observed ones of ellipti¬ 
cal galaxies. In particular, for the most massive, older and 
more metallic, single-burst objects like early-type galaxies, 
our models provide good fits. However, the younger, lower- 
mass counterparts show redder colours that might be at¬ 
tributed to a more extended SFH. These galaxies are most 
likely characterised by the presence of a significant fraction 
of AGB-dominated intermediate-age stellar populations and 
of a large amount of PAH emission corresponding to a rela¬ 
tively recent star formation activity. For most of the galaxies 
showing the reddest ([3.6] — [4.5]) colours, enhanced emission 
in the Spitzer [8] pm band which can be attributed to PAHs 
is indeed observed. Another way to explain these red colours 
would be the contribution of populations of low metallicities, 
for which the Spitzer colour most likely continues to become 
redder and redder. In order to confirm or rule out this op¬ 
tion, we will need more reliable models in this low-metallicity 
regime. 

(iv) We obtained M/L ratios in the K, [3.6]/rm and 
[4.5] pm bands from t he output of our mod elling routines. 
As expected (see, e.g.. iBell fc de Jon^bOOlll . the M/L ra¬ 
tio in all these three bands varies much less depending on 
age and metallicity than when measured in optical filters, 
since light from the RGB stars dominates in the NIR and 
in the IR. Due to the small variations in the (M/Ljs.e ra¬ 
tio with age and meta llicity and thus also with the Spitzer 
([3.6] — [4.5]) colour, iMeidt et al.l ll2014h even suggest a 
parameter-independent M/L ratio of [M/L)z.e = 0.6. This 
value is equal to the mean {M/Ljs.e from our models when 
considering the whole age and met allicity range . This result 
is confirmed by a recent work of iNorris et al.l (l2014li who 
found a variation of the (M/Ljs.e ratio within the same lim¬ 
its measured over a comparable age and metallicity range. 
Whereas the {M/Ljs.e ratio is almost constant as a func¬ 
tion of age and metallicity, it depends significantly on the 
IMF used to calculate the underlying SSP-model. The more 
bottom-heavy the IMF becomes, the more the M/L ratios 
increase. 

(v) Finally, we also compare the behaviour of our SSP 
models to other models available in the literature for our 
specific wavelength range. It is important to note that all of 
them are based on theoretical stellar libraries. The Spitzer 
([3.6] — [4.5]) colours derived from the spectra of our SSP 
models agree alm ost perfectly w i th th e colo ur predictions 
of the models by iMarigo et al.l ll2008ll and iBressan et al.l 
(l2012ll. Also the agreement with the FSPS models of 
IConrov et al.l (120091') is good. The Spitz er colours determined 
from the GALEV (IKotulla et al.ll2009l i models show a large 
offset to the red with respect to the colours measured from 
all other models mentioned so far. This result can be under¬ 
stood by the use of different theoretical stellar atmospheres 
and by the non-consideration of the CO absorption affecting 
the [4.5] ^m band. However, for the Marigo and also for the 
Bressan models, only colour predictions are available, not 
SSP spectra. In the case of the other sets of models, spectra 
can be obtained, but due to their low resolution it is not 
possible to measure spectral features on them. Thus, our 
models are the only ones existing up to now which offer the 
opportunity to obtain line strength indices in the L band. 
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A detailed study of them will be the subject of a following 
paper. 

We thus present the first models in the wavelength range 
between 2.5 and 5 fim which are based on an empirical stel¬ 
lar library. Theoretical stellar spectra were only applied in 
order to fill two gaps in the spectra due to telluric absorp¬ 
tion between « 2.50 and 2.80 /rm and between « 4.20 and 
4.50 /rm, respectively. The very good agreement of the colour 
predictions with comparable theoretical models existing in 
the IR wavelength range between 2.5 and 5/rm as well as 
the possibility to analyse line strength indices renders our 
newly developed models utile. Drawbacks are of course the 
limited range in both age and metallicity due to the limited 
stellar coverage of the IRTF library and the rather average 
fit to particularly the reddest and youngest observed data. 
However, we can state that colours derived from those mod¬ 
els by convolution with the respective filter responses can be 
safely used for all kinds of analyses of observations. Our new 
SSP models will be made available on the MILES webpage 
(http://miles.iac.es). Future work will encompass a detailed 
analysis of line strength indices in the NIR and also to a 
certain extent in the IR as well as c ombining our models 
with the existing MIUSCAT models JVazdekis et al.|[2012l i 
in the optical wavelengths range. 
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APPENDIX A: STELLAR ATMOSPHERIC 
PARAMETERS 

We present our final set of stellar atmospheric parameters 
in the following Table lAll 
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Table Al. Final list of stellar atmospheric parameters: The first column indicates the official names of the stars. Columns 2 and 3 
contain the derived effective temperatures Teff and their errors, respectively. In columns 4 and 5 we present the surface gravities log(p') 
and their errors. Columns 6 and 7 show the metallicities [Fe/H] and their errors. In column 8 we indicate whether a star contributes 
with its full weight (100 per cent) or only with 50 per cent or 25 per cent of it. Colu mn 9 reveals how the effec tive temperatures were 
compiled. An 'A' des ignates that we made use of th e colour-temperature relations of I Alonso et al.l lll996l . [l999l. see Section ES and a 
^CvD^ of the ones of IConrov &; van Dokkuii] ll2012t) . respectively. A means that the respective star was contained in the MILES 
li brary, so that we coul d directly take the value from there. A '1' and a express that we obtained the parameters using the approach 
of ICenarro et af] l|2Q07l. see also Section l5.1|l based on one or seve ral references, respect ively. The abbrevia tions and mark those 
stars for which the respective values were taken from the works of iBergeat et all ll200lf) and IPicklesI lll998h . respectively. 


star name 

Tefl [K] 

error [K] 

log(s) [dex] 

error [dex] 

[Fe/H] [dex] 

error [dex] 

weight 

reference 

HD 31996 

2290 

79 

-0.56 

0.16 

0.26 

0.10 

100 

B 

HD 57160 

3295 

195 

-0.46 

0.16 

0.04 

0.10 

100 

B 

HD 76221 

2645 

24 

-0.46 

0.16 

-0.03 

0.10 

100 

B 

HD 44984 

3035 

10 

-0.46 

0.16 

-0.07 

0.10 

100 

B 

HD 135153 

7049 

47 

1.89 

0.12 

-0.79 

0.10 

100 

1 

HD 13174 

6422 

79 

3.75 

0.10 

-1.90 

0.05 

100 

1 

HD 27397 

7561 

125 

4.05 

- 

-0.91 

0.10 

100 

A 

HD 173638 

7155 

120 

1.81 

0.31 

0.12 

0.05 

100 

+ 

BD +382803 

6597 

- 

1.80 

- 

-1.92 

0.10 

100 

P 

HD 40535 

6805 

- 

3.86 

- 

-0.21 

0.10 

50 

P 

HD 26015 

6855 

47 

4.35 

0.12 

0.17 

0.05 

100 

1 

HD 21770 

6668 

125 

3.72 

- 

-0.65 

0.10 

100 

A 

HD 16232 

6346 

61 

4.54 

0.18 

0.03 

0.09 

100 

M 

HD 87822 

6533 

61 

4.18 

0.18 

0.10 

0.09 

50 

M 

HD 213306 

6051 

- 

1.69 

- 

-0.12 

0.10 

100 

P 

HD 17918 

6701 

125 

3.80 

- 

-0.70 

0.10 

50 

A 

HD 27524 

6622 

61 

4.28 

0.20 

0.13 

0.10 

100 

M 

HD 11443 

6305 

47 

3.67 

0.12 

0.05 

0.10 

100 

1 

HD 126660 

6202 

61 

3.84 

0.18 

-0.27 

0.09 

50 

M 

HD 51956 

6200 

- 

1.77 

- 

-0.09 

0.10 

100 

P 

HD 27383 

6098 

61 

4.28 

0.20 

0.13 

0.09 

100 

M 

HD 114710 

5975 

85 

4.40 

0.28 

0.09 

0.16 

100 

M 

HD 165908 

5928 

85 

4.24 

0.28 

-0.53 

0.16 

100 

M 

HD 185018 

5347 

120 

2.21 

0.31 

-0.24 

0.05 

100 

+ 

HD 20619 

5652 

61 

4.48 

0.18 

-0.28 

0.09 

100 

M 

HD 21018 

5352 

125 

3.07 

- 

0.00 

0.10 

100 

A 

HD 216219 

5727 

85 

3.36 

0.28 

-0.39 

0.16 

100 

M 

HD 10307 

5838 

61 

4.28 

0.18 

0.03 

0.09 

100 

M 

HD 3421 

5383 

120 

2.45 

0.31 

-0.17 

0.05 

100 

+ 

HD 126868 

5651 

120 

3.66 

0.31 

-0.02 

0.05 

100 

+ 

HD 76151 

5692 

61 

4.28 

0.18 

0.08 

0.09 

100 

M 

HD 192713 

5020 

125 

0.05 

0.16 

-0.41 

0.10 

100 

A 

HD 10697 

5610 

120 

3.96 

0.31 

0.11 

0.05 

100 

+ 

HD 165185 

5909 

28 

4.58 

0.07 

0.00 

0.03 

100 

1 

HD 202314 

5007 

92 

1.86 

0.20 

0.13 

0.06 

50 

1 

HD 58367 

4837 

54 

1.94 

0.16 

-0.14 

0.07 

100 

1 

HD 27277 

5081 

125 

2.98 

- 

-0.16 

0.10 

100 

A 

HD 16139 

5102 

- 

2.90 

0.16 

-0.29 

0.10 

100 

P 

HD 25877 

5075 

134 

1.88 

0.18 

0.09 

0.05 

100 

1 

HD 114946 

5171 

61 

3.64 

0.18 

0.13 

0.09 

100 

M 

HD 20618 

5058 

121 

3.10 

0.16 

-0.26 

0.09 

100 

1 

HD 135722 

4847 

85 

2.56 

0.28 

-0.44 

0.16 

100 

M 

HD 75732 

5079 

85 

4.48 

0.28 

0.16 

0.16 

100 

M 

HD 170820 

4604 

61 

1.62 

0.20 

0.17 

0.10 

100 

M 

HD 9852 

4648 

- 

2.82 

- 

0.06 

0.10 

100 

P 

HD 179870 

5246 

125 

1.96 

- 

0.11 

0.10 

50 

A 

HD 124897 

4361 

85 

1.93 

0.28 

-0.53 

0.16 

100 

M 

HD 63302 

4500 

90 

0.20 

0.20 

0.12 

0.09 

100 

M 

HD 10476 

5178 

120 

4.43 

0.31 

-0.16 

0.05 

100 

+ 

HD 23082 

4232 

- 

1.59 

- 

0.04 

0.10 

100 

P 

HD 212466 

3797 

125 

0.01 

- 

0.06 

0.10 

50 

A 

HD 3765 

5051 

120 

4.05 

0.31 

0.03 

0.04 

100 

+ 

HD 187238 

4487 

56 

0.80 

0.10 

0.18 

0.04 

100 

1 
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Table Al - continued 


star name 

T,e [K] 

error [K] 

log(9) [dex] 

error [dex] 

[Fe/H] [dex] 

error [dex] 

weight 

reference 

HD 16068 

4296 

- 

1.80 

0.16 

-0.03 

0.10 

100 

p 

HD 207991 

4152 

- 

2.18 

- 

-0.29 

0.10 

100 

CvD 

HD 181596 

3977 

125 

1.89 

- 

-0.38 

0.10 

100 

A 

HD 36003 

4465 

61 

4.61 

0.18 

0.09 

0.10 

100 

M 

HD 194193 

3785 

125 

1.40 

0.16 

-0.10 

0.10 

100 

A 

HD 209290 

3650 

- 

4.69 

- 

-0.15 

0.10 

100 

CvD 

HD 213893 

3896 

125 

1.62 

- 

-0.04 

0.10 

100 

A 

HD 19305 

3972 

- 

4.67 

- 

0.01 

0.10 

100 

CvD 

BD +60265 

3454 

125 

0.38 

- 

0.18 

0.10 

50 

A 

HD 36395 

3737 

61 

4.90 

0.20 

-1.50 

0.09 

100 

M 

HD 339034 

3278 

125 

-0.25 

- 

-0.07 

0.10 

50 

A 

HD 204724 

3911 

125 

1.45 

- 

-0.71 

0.10 

100 

A 

HD 39801 

3550 

61 

0.05 

0.16 

0.03 

0.10 

50 

M 

HD 219734 

3730 

61 

0.90 

0.18 

0.27 

0.10 

100 

M 

HD 206936 

3746 

121 

-0.04 

0.16 

0.05 

0.09 

50 

1 

HD 10465 

3781 

125 

0.50 

- 

-0.46 

0.10 

100 

A 

HD 95735 

3828 

61 

4.90 

0.20 

-0.20 

0.09 

50 

M 

HD 14488 

3509 

125 

-0.17 

- 

-0.77 

0.10 

50 

A 

HD 28487 

3551 

125 

1.03 

- 

-1.30 

0.10 

100 

A 

HD 14469 

3551 

125 

-0.16 

- 

-0.69 

0.10 

100 

A 

G1 388 

3304 

- 

4.87 

- 

-0.16 

0.10 

100 

P 

G1 268 

3031 

- 

4.76 

0.16 

-0.16 

0.10 

50 

P 

HD 27598 

3647 

125 

1.04 

- 

-0.11 

0.10 

100 

A 

HD 4408 

3571 

125 

1.00 

- 

-0.51 

0.10 

100 

A 

G1 213 

3030 

- 

5.21 

- 

-0.54 

0.10 

100 

CvD 

G1 299 

3111 

- 

5.23 

- 

-0.72 

0.10 

100 

P 

HD 18191 

3250 

61 

0.30 

0.18 

-1.01 

0.10 

100 

M 

HD 196610 

3060 

121 

0.31 

0.16 

-1.02 

0.09 

50 

1 

G1 406 

2698 

- 

4.63 

0.16 

-0.15 

0.10 

25 

P 

HD 207076 

2750 

61 

-0.56 

0.16 

-2.61 

0.10 

100 

M 

HD 113139 

6823 

120 

3.89 

0.31 

0.22 

0.10 

100 

+ 

HD 160365 

6142 

47 

3.10 

0.12 

-0.23 

0.10 

100 

1 

HD 108477 

5172 

125 

3.03 

- 

-0.01 

0.10 

100 

A 

HD 214850 

5437 

- 

4.44 

- 

-0.31 

0.10 

100 

P 

HD 115617 

5531 

61 

4.32 

0.18 

-0.10 

0.09 

100 

M 

HD 122563 

4566 

85 

2.81 

- 

-2.63 

0.16 

100 

M 

HD 101501 

5401 

61 

4.60 

0.18 

-0.13 

0.09 

100 

M 

HD 145675 

5264 

85 

4.66 

0.28 

0.34 

0.16 

100 

M 

HD 2901 

4459 

125 

2.62 

- 

-0.03 

0.10 

100 

A 

HD 132935 

4492 

- 

2.65 

- 

-0.05 

0.10 

100 

CvD 

HD 221246 

4228 

125 

2.32 

- 

-0.24 

0.10 

100 

A 

HD 219134 

4717 

85 

4.50 

0.28 

0.05 

0.16 

100 

M 

HD 216946 

3782 

- 

0.74 

- 

-0.93 

0.10 

100 

P 

HD 237903 

4070 

61 

4.70 

0.18 

-0.18 

0.10 

100 

M 

HD 42581 

3731 

- 

4.66 

- 

-0.16 

0.10 

50 

CvD 

HD 120052 

3664 

125 

1.26 

- 

-0.75 

0.10 

100 

A 

HD 40239 

3349 

121 

0.47 

0.16 

0.03 

0.09 

100 

1 

HD 39045 

3670 

125 

1.12 

- 

-0.99 

0.10 

100 

A 

G1 273 

3108 

- 

5.00 

- 

-0.25 

0.10 

100 

CvD 

HD 94705 

3576 

125 

0.77 

- 

-2.61 

0.10 

100 

A 

HD 108849 

3434 

125 

0.13 

0.16 

-2.27 

0.10 

100 

A 

HD 7927 

7380 

92 

0.80 

0.20 

-0.25 

0.06 

100 

1 

HD 6130 

7041 

120 

2.28 

0.31 

0.06 

0.05 

100 

+ 

HD 89025 

6953 

120 

3.20 

0.31 

0.19 

0.09 

100 

+ 

HD 108519 

6952 

144 

4.16 

- 

-1.51 

0.10 

100 

A 

HD 213135 

6904 

144 

4.17 

- 

-0.53 

0.10 

100 

A 

HD 182835 

7226 

120 

2.12 

0.31 

0.09 

0.05 

100 

+ 

HD 75555 

6458 

47 

3.63 

0.12 

-0.18 

0.10 

100 

1 

HD 186155 

6753 

47 

3.39 

0.12 

-0.25 

0.10 

100 

1 

HD 218804 

6261 

61 

4.05 

0.18 

-0.23 

0.09 

100 

M 

HD 215648 

6167 

61 

4.04 

0.18 

-0.32 

0.09 

100 

M 

HD 201078 

6157 

123 

2.43 

- 

0.13 

0.10 

100 

M 

HD 124850 

6116 

61 

3.83 

0.18 

-0.10 

0.09 

100 

M 
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Table A1 - continued 


star name 

Tea [K] 

error [K] 

log(9) [dex] 

error [dex] 

[Fe/H] [dex] 

error [dex] 

weight 

reference 

HD 102870 

6109 

85 

4.20 

0.28 

0.17 

0.16 

100 

M 

HD 220657 

5878 

47 

3.23 

0.12 

-0.65 

0.10 

100 

1 

HD 111844 

6122 

47 

3.73 

0.12 

-0.74 

0.10 

100 

1 

HD 219623 

6155 

85 

4.17 

0.28 

-0.04 

0.16 

100 

M 

HD 6903 

5684 

47 

2.75 

0.12 

-0.32 

0.10 

100 

1 

HD 176051 

5990 

47 

4.45 

0.12 

-0.18 

0.10 

100 

1 

HD 109358 

5868 

120 

4.28 

0.31 

-0.17 

0.05 

100 

+ 

HD 74395 

5250 

105 

1.30 

0.20 

-0.05 

0.09 

100 

M 

HD 95128 

5813 

61 

4.34 

0.18 

0.04 

0.09 

100 

M 

HD 39949 

5240 

120 

1.23 

0.31 

-0.10 

0.05 

100 

+ 

HD 42454 

5238 

56 

1.10 

0.10 

0.03 

0.04 

100 

1 

HD 219477 

5635 

170 

2.90 

0.16 

-0.12 

0.10 

100 

A 

HD 176123 

5221 

134 

2.29 

0.18 

-0.08 

0.05 

100 

1 

HD 88639 

5280 

125 

3.08 

- 

-0.39 

0.10 

100 

A 

HD 94481 

5302 

54 

2.85 

0.16 

-0.01 

0.08 

100 

1 

HD 190113 

4784 

92 

1.67 

0.20 

0.05 

0.06 

100 

1 

HD 193896 

5050 

170 

2.93 

- 

-0.48 

0.10 

100 

A 

HD 161664 

4405 

125 

1.63 

0.16 

-0.03 

0.10 

100 

A 

HD 333385 

4144 

125 

0.08 

- 

0.27 

0.10 

100 

A 

HD 182694 

5057 

54 

2.94 

0.16 

-0.04 

0.08 

100 

1 

HD 208606 

4701 

92 

1.19 

0.20 

0.12 

0.06 

100 

1 

HD 104979 

4887 

120 

2.55 

0.31 

-0.43 

0.05 

100 

+ 

HD 222093 

4784 

125 

2.85 

- 

0.08 

0.10 

100 

A 

HD 164349 

4446 

61 

1.50 

0.18 

0.39 

0.09 

100 

M 

HD 165782 

4588 

170 

0.54 

- 

0.26 

0.10 

100 

A 

HD 44391 

4670 

120 

0.76 

0.31 

0.12 

0.05 

100 

+ 

HD 100006 

4701 

120 

2.53 

0.31 

-0.17 

0.05 

100 

+ 

HD 91810 

4604 

125 

2.74 

- 

-0.01 

0.10 

100 

A 

HD 36134 

4677 

125 

2.77 

- 

-0.29 

0.10 

100 

A 

HD 25975 

4951 

121 

3.44 

0.16 

0.05 

0.09 

100 

1 

HD 142091 

4788 

120 

3.23 

0.31 

-0.00 

0.05 

100 

+ 

HD 165438 

4862 

97 

3.40 

0.20 

0.02 

0.10 

100 

M 

HD 137759 

4498 

61 

2.38 

0.18 

0.05 

0.09 

100 

M 

HD 114960 

4301 

125 

2.29 

- 

-0.09 

0.10 

100 

A 

HD 99998 

3939 

61 

1.80 

0.18 

-0.31 

0.09 

100 

M 

HD 35620 

4367 

85 

1.75 

0.28 

-0.03 

0.16 

100 

M 

HD 178208 

4473 

170 

2.35 

- 

-0.23 

0.10 

100 

A 

HD 185622 

3848 

- 

0.28 

0.16 

-0.15 

0.10 

100 

P 

HD 201065 

3879 

125 

1.14 

0.16 

-0.49 

0.10 

100 

A 

HD 45977 

4657 

34 

4.27 

0.07 

0.01 

0.05 

100 

1 

HD 120477 

4092 

55 

1.50 

0.16 

-0.56 

0.08 

100 

1 

HD 3346 

3909 

170 

1.19 

0.16 

-0.00 

0.10 

100 

A 

HD 181475 

4024 

79 

1.40 

0.11 

0.00 

0.06 

100 

1 

HD 201092 

3843 

120 

4.50 

0.11 

-0.41 

0.06 

100 

+ 

HD 164136 

6799 

340 

3.96 

- 

-0.30 

0.10 

50 

M 

HD 236697 

3594 

125 

0.52 

- 

0.18 

0.10 

50 

A 

HD 35601 

3617 

120 

0.06 

0.31 

-0.20 

0.05 

100 

+ 

HD 14404 

3640 

125 

0.32 

- 

-0.07 

0.10 

100 

A 

HD 23475 

3863 

170 

0.78 

- 

-0.63 

0.10 

100 

A 

G1 806 

3611 

- 

4.71 

- 

-0.19 

0.10 

100 

GvD 

HD 204585 

3428 

- 

1.00 

- 

-1.46 

0.10 

100 

GvD 

HD 19058 

3544 

125 

0.99 

- 

-0.96 

0.10 

100 

A 

HD 214665 

3397 

125 

0.91 

- 

-1.40 

0.10 

100 

A 

HD 156014 

3236 

- 

0.73 

0.16 

-2.61 

0.10 

100 

P 

HD 175865 

3420 

61 

0.50 

0.18 

0.14 

0.10 

100 

M 

G1 51 

2794 

- 

5.36 

- 

-0.15 

0.10 

100 

GvD 

G1 866 

2951 

- 

5.27 

- 

-0.15 

0.10 

100 

P 

HD 190323 

5715 

56 

0.40 

0.10 

-0.24 

0.04 

100 

1 

HD 18474 

4899 

120 

2.36 

0.31 

-0.30 

0.05 

100 

+ 

HD 92055 

2965 

57 

-0.56 

0.16 

-0.06 

0.10 

100 

B 

HD 76846 

4687 

125 

2.10 

0.16 

-0.11 

0.10 

100 

A 

G1 381 

3319 

- 

4.64 

- 

-0.00 

- 

100 

GvD 

G1 581 

3197 

- 

4.57 

- 

-0.00 

- 

100 

GvD 

HD 179821 

4706 

- 

0.34 

- 

-0.00 

- 

100 

A 



